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I„    INTRODUCTION 

In  recent  years,  increasing  attention  has  been  given  to  the  problem  of 
homogeneous  nucleation  and  condensation  of  supersaturated  vapors. 

Particular  interest  has  arisen  in  connection  with  high  temperature, 
low  vapor  pressure  fluids,  capable  of  high  Carnot  efficiencies  in  power 
plant  application.  The  possible  space  application  of  systems  employing 

alkali  metals,  and  recent  reports  of  experiments  involving  the  flow  of 

1+3  5I1 

sodium  and  potassium  through  nozzles   '   underlines  the  need  to  understand 

the  condensation  behavior  of  these  gases  well  enough  to  provide  a  basis  for 
optimum  turbine  and  cycle  design. 

Early  investigations  in  cloud  chambers  showed  that,  after  initial 
expansions  and  settling  had  removed  dust  particles  and  other  nuclei,  sub- 
sequent expansions  into  the  saturated  region  would  result  in  the  spontaneous 
formation  of  a  condensation  cloud  at  some  limiting  expansion  ratio.   In  the 
absence  of  other  nucleation  sites,  it  was  concluded  that  such  clouds  were 
the  result  of  homogeneous  nucleation  of  the  super-saturated  vapor.   Similar 
condensation  occurs  in  vapors  expanded  in  nozzles,  where  rates  of  expansion 
are  generally  much  more  rapid  than  those  of  cloud  chambers. 

Despite  extensive  literature  on  the  subject,  there  exists  little  agree- 
ment as  to  the  correct  formulation  of  nucleation  and  drop  growth  rates,  and 
within  the  theoretical  framework  there  exist  many  other  uncertainties  regard- 
ing the  fluid  properties  pertinent  to  numerical  evaluation.   Experimental 
work  is  relatively  limited  in  terms  of  range  of  test  conditions,  and  numbers 
of  fluids  tested.   Although  a  variety  of  gases  have  been  studied  in  cloud 
chamber  expansions,  applicable  and  informative  studies  of  nozzle  expansions 
are  somewhat  sparse  outside  of  steam  and  airborne  water  vapor  for  which  both 
cloud  chamber  and  nozzle  data  is  extensive.  *   '   '        This  data  as  well 
as  the  cloud  chamber  data  for  several  other  vapors  tends  to  substantiate 
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so-called  "classical"  '   '   nucleation  theory,  as  opposed  to  revised  models 
of  nucleation  theory  offered  by  more  recent  authors,  as  a  better  means  of 
predicting  condensation  behavior „   However,  there  is  doubt  regarding  the 
extent  to  which  cloud  chamber  data  may  be  applied  to  the  determination  of 
nucleation  rates,  and  too  little  data  outside  of  water  vapor  has  been  obtained 
for  flow  in  nozzles  to  allow  a  comprehensive  test  of  theory  over  a  variety 
of  gases o 

The  purpose  of  this  work  has  been  to  extend  the  limits  of  existing 
experimental  data  involving  condensation  in  nozzles,  through  experimentation 
with  C0_,  to  study  by  means  of  a  computational  program  the  predictions  of 
existing  theory  and  make  some  judgments  regarding  the  importance  of  the 
various  physical  parameters  and  uncertainties  pertaining  thereto,  and  to 
reduce  these  uncertainties  to  a  point  where  a  valid  comparison  of  existing 
theory  and  experimental  results  can  be  made  for  CO  „ 
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II.   NOZZLE  FLOW  PHENOMENA 

Converging-diverging  nozzles  lend  themselves  particularly  nicely  to 
the  study  of  condensation  "behavior  of  gases.   Following  are  some  of  the 
aspects  of  gas  behavior  in  nozzle  flow  and  features  which  recommend  it  as  a 
means  of  experimental  study  of  condensation. 

A„   Isentropic  Effects 

It  is  general  knowledge  that  for  sufficiently  high  pressure  ratios,  an 
ideal  gas  will  expand  isentropically  in  a  converging-diverging  nozzle, 
becoming  sonic  at  the  point  of  minimum  flow  area  and  supersonic  downstream, 
with  continuously  increasing  velocities  and  decreasing  pressures,  densities 
and  temperatures  all  along  the  expansion.   A  non-ideal  gas  behaves  similarly, 
and  depending  upon  the  stagnation  conditions,  the  extent  of  the  expansion, 
and  the  properties  of  the  vapor  being  tested,  the  isentrope  may  approach 
saturation  conditions  of  pressure  and  temperature.   As  illustrated  on  the 
T-S  diagram  of  Figure  1,  it  is  possible  for  conditions  to  exist  for  certain 
gases  when  the  opposite  may  be  true  and  an  isentropic  expansion  will  result 
only  in  increasing  the  degree  of  superheat  of  the  vapor  (i.e.  as  with  many  of 
the  FREONS).   However,  in  most  cases  it  is  a  simple  matter  to  devise  stagna- 
tion conditions  which  result  in  the  rapid  arrival  of  the  gas  to  conditions 

of  super saturation  during  its  isentropic  expansion.   For  a  supersonic  nozzle 

-1* 

a  few  inches  in  length  expansion  time  is  of  the  order  10   seconds. 

B.   Wall  and  Boundary  Layer  Effects 

In  all  real  gases  there  are  departures  from  ideal  non-viscous  behavior 
which  must  be  considered.   Fortunately,  in  nozzle  flow  these  effects  are 
easily  accounted  for,  provided  fully  developed  flow  does  not  occur.   Viscous 
effects  are  accounted  for  by  considering  the  nozzle  flow  to  be  an  isentropic 
expansion  contained  between  wall  boundary  layers.  The  boundary  layers  are 
greatly  influenced  by  the  shape  of  the  pressure  gradient  and  are  generally 


depressed  to  a  very  small  mimimum  value  near  the  most  critical  point,  the 
geometric  throat 0   Further 9   viscous  and  recovery  effects  on  the  nozzle  walls 
will  keep  the  wall  temperature  close  to  stagnation  value  in  the  absence  of 
very  strong  wall  cooling  or  heating,  and,  at  decreased  downstream  pressures , 
is  adequate  guarantee  that  there  will  be  no  condensation  of  liquid  on  the 
walls  which  could  interfere  with  the  mid-stream  process „ 

In  experimentation  where  nozzle  size  might  be  limited  by  mass  flow  or 
other  restrictions,  considerations  of  Reynold tfs  Numbers  and  boundary  layer 
thicknesses  which  may  be  expected  in  the  throat  region  are  important  since 
large  boundary  layers  in  relation  to  nozzle  size  could  invalidate  the 
assumption  of  an  isentropic  region  of  flow0  A  means  of  estimating  the 
boundary  layer  thickness  in  the  throat  region  of  a  supersonic  nozzle  was 
devised  (Appendix  E)  and  showed  good  agreement  with  interferometric  photos 
that  the  boundary  layers  involved  in  the  CO  experiments  of  this  investiga- 
tion were  about  0OO611  at  the  nozzle  throat  0   In  all  nozzles  tested  this  was 
less  than  11$  of  the  throat  width „  Additional  discussion  of  boundary  layer 
treatment  in  this  investigation  is  included  later0 

The  isentropic  region  between  boundary  layers  may  be  described  by  one- 
dimensional  gas  dynamics a  and,  by  means  of  measurement  of  some  vapor 
property  (such  as  static  pressure)  along  the  nozzle,  a  complete  time- 
history  of  the  expansion  may  be  obtainedo 

Condensation  in  a  nozzle  expansion  has  been  shown  by  many  investiga- 
tions to  have  a  large  influence  on  the  thermodynamic  properties  of  the 
fluid o   In  particular,  the  static  pressure  records  large  variation  from  the 
non»condensing  profile  as  a  result  of  the  heat  released  by  moisture  forma- 
tion o   For  the  addition  of  heat  to  an  ideal  gas,  it  may  be  shown  that 

§L     -  YM2    rdA    dQ_i 
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For  water  vapor,  it  is  easily  estimated  that  at  M  =  2,  a  formation  of 
1%   moisture  at  constant  area  is  accompanied  by  about  1%   variation  in  static 
pressure,,  As  is  later  shown,  it  is  possible  to  detect  the  onset  of  condensa- 
tion at  about  al%   presence  of  moisture  through  simple  static  pressure 
measurements o 

It  may  be  seen  from  the  above  equations  that  pressure  variation  is 
increased  by  establishing  conditions  such  that  condensation  occurs  with 

Mach  Number  close  to  unity,, 

6k  3 

Willmarth  and  Nagamatsu,   and  Arthur,  have  shown  that  the  ratio  P/P' 

o 

plotted  against  the  streamwise  variation  of  P'/P  is  an  even  more  sensitive 
indicator  of  condensation,  where  P  is  the  measured  static  pressure,  P*  is 
the  measured  total  pressure  behind  a  normal  shock,  and  P  is  the  original 
stagnation  pressure , 

An  important  consideration  in  any  system  designed  to  test  homogeneous 
nucleation  theory  is  the  possible  effects  of  contamination^,   As  will  be  shown 
later,  nozzle  expansions  from  superheated  to  supersaturated  state  proceed  in 
such  an  extremely  short  time  that  spontaneous  nucleation  will  produce  nuclei 
many  orders  of  magnitude  greater  in  number  than  those  which  could  be  present 
in  the  form  of  dust  or  other  particles  and  the  amount  of  moisture  formed  on 
such  contaminants  is  entirely  negligible „  As  will  also  be  mentioned  later, 
an  exception  to  this  is  the  case  of  a  pre-condensing  vapor „ 

The  proceeding  observations  indicate  that,  with  fairly  simple  instru- 
mentation of  a  nozzle  system,  it  is  possible  to  obtain  a  complete  history  of 
a  gas  expansion  and  any  condensation  occurring  within,  with  the  reasonable 
assumption  that  no  outside  influences  exist  which  would  invalidate  testing 
of  homogeneous  nucleation  and  drop  growth  theory  in  the  same  system, 

Co   Test  Vapor  Considerations 

The  selection  of  a  vapor  for  testing  is  subject  to  the  following 


considerations „   First,  a  large  number  of  vapors  are  toxic,  reactive  to 
instrumentation  materials,  or  do  not  exist  as  vapors  at  conditions  most 
convenient  for  laboratory  testing.,   Second,  the  state  relations  and  other 
properties  of  many  are  so  poorly  defined,  that  application  to  any  meaningful 
test  of  theory  would  be  ill-foundeda   Neither  of  these  considerations  would 
be  restrictive  if  the  condensation  data  for  a  particular  vapor  were  parti- 
cularly needed,  or  if  a  test  of  theory  were  incidental  to  the  experimental 
data»   Third,  two  particular  properties  are  desirable  in  a  test  vapor,  both 
of  which  are  more  common  to  gases  of  low  molecular  weight .  The  first  has 
been  mentioned  previously 9  and  involves  the  shape  of  the  coexistence  line  in 
relation  to  the  isentrope„  A  preferred  test  gas  has  a  coexistence  line 
having  rapid  convergence  with  the  isentrope„   Also,  the  measure  of  heat 
release  resulting  from  condensation  compared  with  that  needed  to  alter  its 
vapor  properties  is  important „   This  may  be  expressed  through  the  parameter 
h  /C  ,  for  which  typical  values  for  various  gases  are  presented  in  Table  1„ 
Variations  in  the  gas  properties  caused  by  condensation  are  increased  with 
larger  values  of  this  parameter „  A  more  applicable  use  of  these  values  is, 

as  seen  from  Equation  1,  the  form  h„  /C  T  at  the  test  conditions  of  interest, 

fg  P 

For  vapors  such  as  alkali  metals,  higher  values  of  h  /C  are  decreased  in 

1 8  P 

significance  due  to  the  high  temperatures  of  the  gas» 

C0_  was  selected  for  this  study  because  it  satisfies  all  the  foregoing 
criteria  and  in  addition  may  be  obtained  in  high  purity  at  relative  low 
cost j,  thus  allowing  its  use  in  a  system  not  necessitating  recovery  after 
expansion o   Also,  such  physical  properties  as  CO  vapor  pressure,  liquid 
surface  tension,  and  solid  and  liquid  density,  necessary  in  the  application 
of  theory,  are  well  established  and  readily  available  in  numerous  publica- 
tions of  physical  property  datac 


Discussion  of  the  application  of  these  properties,  and  details  of  the 
theoretical  mechanism  and  behavior  of  condensation  are  presented  later „ 

The  question  could  be  raised  concerning  possible  non-equilibrium  of 

30 
energy  modes  of  CO  during  such  a  rapid  expansion „   Kantrowitz   has  shown 

that  this  would  not  be  a  problem  with  CO  in  this  case  since  1) vibrational 
energy  modes  are  not  activated  at  the  low  temperatures  of  experimental  test 
conditions  and  2)the  relaxation  time  constant  of  vibration  is  considerably 
less  (=10°"  sec)  than  that  of  the  gas  expansion.,   Camac   has  recently 
formulated  the  pressure-temperature  dependence  of  CO  vibrational  relaxa- 
tion time  and  his  expression  also  yields  times  of  the  order  10   sec„  in 
the  region  of  experimental  testings 
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III.  EXPERIMENTAL  INVESTIGATIONS 

A.   Test  Apparatus 

Industrial  Welding  CO  having  very  high  (moisture  content  less  than 
.00*45/0  by  weight)  purity  was  introduced  via  a  manifold  from  a  bank  of  nine 
50  lb  cylinders  to  a  needle  valve  used  for  regulation  of  flow  rate  (See 
Figure  2).   Thence  the  flow  entered  a  stagnation  tank  and  then  the  test 
nozzle.   Exhaust  was  to  atmosphere  in  earlier  tests  and  later  to  a  steam 
ejector  system  maintaining  about  23"  Hg„  of  vacuum.   This  modification 
facilitated  lower  minimum  stagnation  pressures  without  the  interference  of 
compression  shock  systems  in  the  supersonic  test  section.   Nine  cylinders 
were  employed  in  order  to  obtain  a  large  mass  flow  source  providing  nearly 
constant  conditions  to  the  nozzles  for  extended  periods  of  flow,  despite 
the  quasi-steady  nature  of  the  "blow-down"  method. 

The  stagnation  tank  contained  a  calibrated  pressure  gauge  and,  extend- 
ing into  the  flow  just  upstream  of  the  nozzle  entrance,  a  copper-constantan 
thermocouple  using  an  ice  reference.   A  100  inch  25-tube  mercury  filled 
manometer  board  was  used  in  connection  with  static  wall  taps  to  determine 
pressure  readings  along  the  nozzle.   Pressure  data  was  recorded  by 
polaroid  photograph  of  the  manometer  board.   Test  time  durations  were  of  the 
order  of  1  minute. 

The  stagnation  tank  had  a  removable  end  plate  and  contained  screens 
which  allowed  the  packing  of  large  amounts  of  steel  wool  inside  a  nylon 
filter  bag  upstream  of  the  nozzle  entrance.  (Figure  5)   When  high  stagnation 
temperatures  were  desired  (i.e.  non-condensing  tests)  an  electric  heating 
tape  wrapped  around  the  stagnation  tank  was  energized  and  the  steel  wool 
used  for  heat  storage  and  thermal  inertia  to  decrease  the  rate  of  the 
temperature  transient.   When  low  temperatures  were  desired,  the  steel  wool 
was  removed  and  the  heating  tape  not  used. 
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B.    Test  Nozzles 

Two  test  nozzles  were  employed  for  experimental  measurements, 
henceforth  referred  to  as  Nozzle  I  and  Nozzle  II . 

Nozzle  I  appears  schematically  in  Figure  7  and  in  photos  in  Figures  6 
and  9.   It  featured  two  removable  parallel  side  walls  which  were  plates 
of  optical  glass,  used  for  interferometric  measurements.   By  exchanging 
one  of  the  plate,  for  an  identical  metal  one  containing  pressure  taps, 
static  pressure  profiles  could  be  obtained  and  the  flow  observed  simulta- 
neously through  the  remaining  glass  wall.   Slotted  "O-rings"  provided 
sealing  of  the  side  walls.   The  nozzle  was  made  wide  in  comparison  to  its 
height  in  order  to  minimize  any  contribution  to  error  caused  by  boundary 
layer  effects  of  the  side  walls  on  the  interferometric  measurements.   This 
error  was  later  estimated  to  be  less  than  .2%   of  the  measured  densities. 

The  nozzle  shape  was  also  intended  to  provide  nearly  2-dimensional  flow, 
and  in  fact,  interferometric  photos  to  be  discussed  later  show  the  flow  to 
be  effectively  1-dimensional  between  boundary  layers.   A  region  of  isentropic 
flow  was  thus  obtained  between  the  boundary  layers,  which  were  of  thickness 
=  „005"  -  .007"  at  the  throat  (Appendix  E). 

Pressure  Taps  in  Nozzle  I  were  spaced  .2"  apart,  and  later  during  the 
program,  additional  taps ,  resulting  in  .1"  spacing,  were  added  in  the  test 
section  region. 

Nozzle  II  is  illustrated  in  Figures  8  and  10  and  featured  hinged 
upper  and  lower  nozzle  walls  which  could  be  placed  at  arbitrary  angles 
between  the  fixed  parallel  side  walls.   This  facilitated  variable  expan- 
sion rates  in  the  test  section,  depending  on  the  divergence  angle  of  the 
test  section,  and  it  was  hoped  to  thus  establish  a  "best"  angle  of 
divergence  to  obtain  the  most  informative  pressure  profiles.   Since 
pressure  variations  due  to  heat  addition  approach  maximum  when  Mach 
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Number  approaches  1.00,  it  was  expected  that  this  result  would  be  obtained 
for  a  very  low  angle  of  divergence . 

Static  pressure  taus  located  in  the  side  wall  were  located  .25  inches 
apart  for  early  runs,  after  which  a  new  test  plate  with  taps  .1  inches 
apart  in  the  test  region  was  substituted. 

C .   Experimental  Procedures 

1.    Pressure-Density  Measurements 

Tests  were  conducted  using  two  persons  and  the  following  procedure. 

1)  The  needle  valve  was  opened  until  the  approximate  value  of 
desired  stagnation  pressure  was  obtained. 

2)  One  observer  maintained  a  null-balance  on  the  thermocouple- 
potentiometer  circuit,  while  the  other,  manning  a  tripod-mounted  polaroid 
camera,  observed  the  manometer  board. 

3)  When  the  manometer  board  reached  a  steady  state,  a  polaroid 
photo  was  taken  of  the  manometer  board  by  the  second  observer  while  the 
first  recorded  the  stagnation  gauge  pressure  and  terminated  the  null- 
balancing  procedure  on  the  thermocouple  potentiometer. 

h)    The  reading  of  the  potentiometer  at  time  of  photograph  was 
converted  to  a  stagnation  temperature  reading,  and  the  polaroid  display  of 
the  pressure  data  developed,  and  all  pertinent  data  involving  the  test 
run  affixed  to  it  (including  local  barometric  pressure). 

Interferometer  runs  were  made  in  a  similar  manner,  except  for  the  use 
of  5"  x  7"  sheet  film  in  the  built-in  camera  of  the  interferometer,  and 
requiring  only  the  potentiometer  null-balance  prior  to  recording  data. 

Discussion  of  interferometry  theory  and  the  methods  applied  in  this 
investigation  plus  related  references  are  presented  in  Appendix  G. 
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2 .   Stagnation  Pressure  Measurements 

The  stagnation  pressure  behind  a  normal  shock  was  measured  by  means 
of  a  .030"  O.D.  pitot  tube  extending  upstream  from  a  micrometer  mounted  on 
the  elbow  downstream  of  the  nozzle. (See  Figure  6)   All  measurements  were 
made  on  Nozzle  I,  by  visually  aligning  (with  the  aid  of  lighting,  and 
machinist's  templates)  the  tip  of  the  tube  with  the  center  of  the  corres- 
ponding wall  tap  with  which  static  pressure  was  measured.   Since  the 
location  of  the  shock  was  actually  located  a  small  distance  upstream,  a 
correction  was  made  to  the  static  pressure  reading  corresponding  to  this 
distance,  (about  .63  tube  diameters  at  Mach  I.65  and  about  .30  tube  dia- 
meters at  Mach  2.U5,  decreasing  further  as  Mach  Number  increases.) 
Upstream  stagnation  pressures  were  measured  with  the  calibrated  gauge  in 
the  stagnation  tank.   As  in  the  case  of  the  static  pressure  profile  measure- 
ments, photographs  of  the  manometer  board  were  taken  simultaneously  with 
stagnation  temperature  and  pressure  measurements,  after  steady  pressure 
readings  and  a  null  potentiometer  reading  were  obtained. 

Total  pressure  measurements  were  made  with  all  manometer  tubes  dis- 
connected except  the  following  l)pitot  tube  2)static  wall  tap  at  tip  of 
pitot  tube  3) throat  static  wall  tap  (recorded  only  for  supplementary 
purposes ) . 

Discussion  and  estimates  of  possible  experimental  error  are  contained 
in  Appendix  I . 

D .   Non-Condensing  Measurements 
1 .   Real  Gas  Behavior 

It  became  necessary  to  consider  the  departure  from  perfect  gas  behavior 
of  CO  near  the  very  outset  of  this  investigation,  since  one  of  the  objects 
of  experimental  work  was  to  compare  the  agreement  of  measured  density  pro- 
files with  measured  pressure  profiles.   Failure  to  consider  the  proper 
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p_v-T  relation  for  reduction  of  interferometric  data  resulted  in  consider- 
able error  in  the  density  profile.   It  also  introduced  doubt  as  to  how  the 
proper  density  ratios  should  be  converted  to  the  corresponding  isentropic 
pressure  ratios  for  purposes  of  making  comparison  between  measurements. 
Existing  tabulated  properties  of  CCL  also  showed  that  the  degree  of 

departure  from  perfect  gas  behavior  for  CO   (as  reflected  by  the  compressi- 

PvN 
bility  factor,  z  =  — )  increased  markedly  near  the  saturated  region,  and 

that  the  values  of  compressibility  factor  could  easily  be  as  low  as  .9 
during  an  expansion.   Also,  existing  tables  did  not  extend  to  pressures  and 
te/mperatures  low  enough  to  make  adequate  determination  of  what  the  real 
gas  behavior  of  CO  is . 

Since  a  meaningful  test  of  the  nucleation  and  drop  growth  predictions 
which  were  to  follow  must  rest  on  the  most  precisely  accurate  information 
available  regarding  the  gas  dynamics,  and  since,  as  will  be  shown,  the  pre- 
dictions are  quite  sensitive  to  errors  in  this  respect,  it  was  decided  to 
establish  the  real  gas  behavior  of  CO  through  use  of  an  applicable  equa- 
tion of  state. 

Unfortunately,  though  the  literature  abounds  with  work,  much  of  it 

i_.  ,  *  „•  U,11,2U,26,28, 

recent,  on  high  pressure  and  temperature  properties  of  CO  , 

33  51     &  59 

it  was  possible   to   find   only  one   equation   of  state  that    is 

applicable   in   the  lower  P-T  range   of  these   experiments,   that  being 

Plank's      (1929),   which   fits   the  data  of  Lowry   and  Erikson,        and  Maass   and 

Mennie        very  well.      Comparison   of  Plank's   Equation  of  State  and  of  their 

Pv 
data  is    shown  in  Figure   In,    showing  maximum  disagreement   of  about    .1%   of  — , 

Kl 

For  comparison,  the  equation  of  State  of  Martin  &  Hou,  and  Benedict,  Webb, 
and  Rubin, applied  well  outside  their  applicable  region  to  the  area  of 
experimental  interest,  vary  by  as  much  as  10$  from  Plank's. 
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Plank's   equation   is  statedand  was   applied  as    follows: 

RT        (.0825  +    .001265Tp)    x  103  ,    . 

v  =  — -  +  = * (2) 

P  (.OIT)^ 

Using  the  relations: 

dp    S  p 

ah  ■  (H'PdT  +  (!pVp  <">* 

■  vT  +  ^  Vp  (Mb 

1  3^ 

where  t  =  — ,  it  is  possible  to  use  the  zero  pressure  value   for  C   in  the 

result         T         p  ./   \ 

h  =  /   C  dT  +  /  [4^  ]  dp  ,M 

J0   P       JOl  dT  Jp  *  (U)c 

The  right  hand  term  may  be  evaluated  from  the  equation  of  state  and 
hence  it  is  possible  to  determine  enthalpy  at  any  desired  pressure  and 
temperature.   To  determine  an  isentrope,  equation  (3)  is  used  and  the  change 
in  h  calculated  for  an  arbitrary  change  in  p.   Thence,  knowing  new  values 
of  p  and  h,  the  corresponding  value  of  T  may  be  determined,  either  through 
the  use  of  pre-calculated  enthalpy  tables,  or  by  iterative  integration  of 
equation  (h)    .   It  is  then  a  simple  matter  to  calculate  the  new  value  of 
density  from  equation  (2)  and  to  thus  determine  a  series  of  pressure, 
temperature,  and  density  ratios  along  an  isentrope  for  any  given  set  of 
stagnation  conditions . 

Table  2  contains  values  of  enthalpies  calculated  in  this  manner. 
Agreement  within  .1°  to  .3°K  occurs  at  the  juncture  of  those  presented 
by  Din.15 

In  addition  to  the  comparison  of  pressure  and  density  measurements,  it 
was  decided  to  apply  an  additional  experimental  test  to  the  ability  of  the 
Plank  Equation  of  State  to  predict  the  behavior  of  CO  gas  dynamics  in 
regions  of  low  pressure  and  temperature  yet  below  its  known  region  of 
applicability,  (0°C  to  -T0°C) .   Measurement  of  the  normal  shock  pressure 
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characteristics  were  made  and  compared  with  the  predictions  of  "real"  gas 
dynamics.   For  this  purpose,  it  was  necessary  to  include  a  computation 
scheme  for  solution  of  the  continuity,  momentum,  energy,  and  state  equa- 
tions across  a  normal  shock  and  to  return  along  a  new  isentrope  to  the  new 
stagnation  conditions .   A  listing  of  this  computation  procedure  is  given  in 
Appendix  H. 

2.   Pressure-Density  Measurements 

Photographs  of  two  interferometer  fringe  patterns  for  non-condensing 
flow  are  seen  in  Figures  11  and  12  and  the  resulting  densities  given  in 
Table  3  and  k.      The  photographs  are  marked  by  an  extended  region  of  near 
verticle  interference  fringes  between  the  upper  and  lower  nozzle  walls, 
indicating  that  this  region  is  in  1-dimensional  flow,  and  that  the  boundary 
layers  are  quite  small  (approximately  .006"  at  the  throat). 

Results  of  the  two  sets  of  measured  density  and  the  corresponding 
pressure  profiles  of  Nozzle  I  are  shown  in  Figures  15a,  15b  and  l6a,  l6b. 
Each  comparison  has  been  illustrated  in  two  ways.   Figures  15a  and  l6a  show 
the  density  and  pressure  profiles  in  relation  to  nozzle  position.   Pres- 
sures have  been  directly  measured  and  plotted  at  their  tap  locations. 
Measured  density  ratios  have  been  converted  to  their  corresponding  pressure 
ratios  through  use  of  the  real  gas  isentropic  values  based  on  the  experimen- 
tal stagnation  conditions. 

Figures  15b  and  l6b  show  the  same  data  plotted  in  somewhat  different 

form,  illustrating  the  theoretical  isentropic  values  of  P/P  vs.  p/p   and 

those  values  of  P/P  and  p/p   which  were  measured  at  identical  points  in 

o        o 

the  nozzle.   On  Figure  l6b,  the  corresponding  perfect  gas  conversion  from 
density  to  pressure  (based  on  zero  pressure  value  of  C  /C  )  is  also  shown. 

\J  V 
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Figures  15  and  16  show  very  good  agreement  between  measured  pressures 
and  densities,  and  those  predicted  from  the  Plank  Equation  of  State.   The 
fact  that  the  conversion  of  density  to  pressure  is  not  very  significantly 
affected  by  the  use  of  the  equation  of  state  rather  than  perfect  gas 
relations,  is  indicated  by  the  close  proximity  of  the  perfect  gas  line  in 
Figure  l6b,  where  the  applied  ratio  of  specific  heats  is  the  zero  pressure 
value  at  stagnation  temperature,  y  =  1«28.  However,  it  should  be  emphasized 
that  even  this  use  of  the  perfect  gas  approximations  is  preceeded  by  the 
use  of  real  gas  p-v-T  relations  in  the  reduction  of  interferometric  fringe 
shift  data  to  obtain  values  of  density.   The  only  noticeable  departure  from 
exact  real  gas  agreement  consistent  to  both  tests  occurs  near  the  throat 
region  where  the  measured  pressures  are  slightly  higher  {-   1.8"  Hg)  than 
those  values  based  on  measured  density.   This  departure  is  concluded  to  be 
due  to  departure  from  1-dimensional  flow  over  the  boundary  layer  of  the 
lateral  side  walls  in  this  region  and  an  approximate  numerical  estimate  of 
this  effect  in  good  agreement  with  the  noted  error  is  made  in  Appendix  F. 

Attempts  to  apply  the  interferometer  to  direct  study  of  condensation 
were  frustrated  by  increasing  fluctuation  and  blurring  of  the  interference 
fringes  as  stagnation  pressure  was  increased.   This  effect  was  apparently 
due  to  thermal  turbulence.   Above  pressures  of  85  psia,  complete  blurring 
resulted  upstream  of  the  throat,  although  the  interference  bands  again 
became  visible  downstream,  where  acceleration  of  the  flow  dampens  the 
effects  of  turbulence.   Consequently,  no  density  profiles  involving  conden- 
sation were  obtained  although  photographs  of  the  condensation  ice  cloud 
showing  the  interference  fringes  within  were  obtained. 
3.   Normal  Shock  Total  Pressure  Measurements 

A  more  sensitive  test  of  the  Plank  Equation  of  State  predictions 
versus  those  of  perfect  gas  was  made  by  applying  each  to  the  gas  dynamics 
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of  normal  shock  behavior  and  comparing  experimental  results  with  theoretical 
predictions  based  on  the  experimental  stagnation  conditions  of  each  test. 
Results  of  experimental  normal  shock  measurements  are  shown  in 

Figure  IT,  using  the  display  of  Willmarth  and  Nagamatsu,  and  Arthur,  who 

6U  3 
made  similar  measurements  for  Nitrogen  and  Helium.   '   The  predictions 

computed  with  Plank's  Equation  of  State  are  shown  in  Figure  18 .   Test 
aeries  No.  2  was  taken  with  a  great  deal  of  regard  for  better  accuracy  in 
pitot  tube  positioning  and  establishment  of  steady  state  readings  prior  to 
taking  data,  with  a  resulting  decrease  in  the  amount  of  scatter.   One- 
dimensional  perfect  gas  predictions  are  shown  for  reference  in  the  lines  of 
constant  y»  the  appropriate  zero  pressure  value  for  most  of  the  test  con- 
ditions being  about  y   =  1.28.   It  is  seen  that  the  real  gas  predictions  are 
considerably  more  accurate  than  those  of  perfect  gas  in  predicting  the 
normal  shock  characteristics  of  CO  ,  although  there  still  exists  some 
departure  of  the  experimental  results  from  those  predicted  with  the  Plank 
Equation  of  State.   The  experimental  results  indicate  behavior  very  close 
to  that  of  a  perfect  gas  having  a  value  of  y   =  l.UO.   However,  this  is 
misleading,  since  treatment  of  CO  as  a  perfect  gas,  using  this  value  of 
ratio  of  specific  heats,  does  not  successfully  correlate  the  pressure- 
density  measurements,  or  correctly  predict  local  temperatures  along  an 
isentrope.   In  terms  of  the  measured  properties,  the  experiment  indicated 
about  2%   larger  values  of  stagnation  pressure  behind  the  normal  shock  than 
those  predicted. 

Despite  the  above  mentioned  lack  of  complete  agreement  between  Plank's 
Equation  of  State  and  measured  normal  shock  pressure  characteristics,  it 
appears  that  such  a  "real  gas"  treatment  of  the  gas  dynamics  is  a  con- 
siderably better  description  than  perfect  gas  and  so  warranted  application 
to  the  prediction  of  local  conditions  along  the  isentrope. 
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Influence  of  Equation  of  State  Upon  Isentropes 

Hie  influence  of  real  gas  treatment  to  the  computation  of  CO  isentropes 
has  quite  a  significant  effect  upon  the  values  of  local  temperature  calcu- 
lated.  Figures  19  and  20  compare  two  typical  perfect  gas  isentropes  with 
those  predicted  using  Plank's  Equation  of  State.   These  curves  indicate  that 
an  error  of  5°K  could  easily  occur  by  attempting  to  describe  the  i sent rope 
by  a  constant  ratio  of  specific  heats,  and  that  a  decrease  in  stagnation 
temperature  could  cause  errors  of  possibly  10°K  if  one  isentrope  were  applied 
to  all  experimental  test  conditions. 

The  resulting  uncertainty  in  predicted  nucleation  rate,  depending  upon 
location  within  the  supersaturated  region,  can  be  seen  by  referring  to  the 
classical  nucleation  rate  profiles  for  CO  illustrated  in  Figure  21. 
Errors  of  15  orders  of  magnitude  could  easily  occur. 

E.    Condensation  Measurements 

1.   Pressure  Profiles  and  Local  Conditions  at  Onset  of 
Condensation 

Results  of  experimental  measurements  involving  condensation  in 
Nozzle  I,  pressure  tap  plates  1  &  2  are  shown  in  Figures  22  and  23. 
Results  of  Nozzle  II,  pressure  tap  plate  1  are  shown  in  Figure  2U,  and  of 
the  same  nozzle  using  pressure  tap  plate  2  in  Figures  25  through  30. 
Figures  2k,    25,  26,  27  and  28  are  for  Nozzle  II  at  an  included  divergence 
angle  of  approximately  2.2°  while  Figures  2Q  &  30  are  for  Nozzle  II  at  an 
angle  of  0.80°.   For  convenience  and  due  to  manometer  board  limitations, 
pressure  ratios  are  expressed  in  terms  of  the  pressure  at  the  tap  nearest 
the  throat,  but  may  be  easily  converted  to  total  pressure  ratio.   Data 
for  each  Figure  is  adjoined  in  Tables  5  through  13. 

All   experimental  profiles    showed  a  marked  occurrence  of  pressure 
variation   from  non-condensing  values   at   the  point   of  occurrence   of 
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condensation  (verified  visually  in  Nozzle  I).   By  noting  the  pressure  ratio 
at  which  condensation  occurs,  it  is  possible  to  determine  the  local  condi- 
tions at  occurrence  using  the  isentropic  relations  for  the  equation  of 
state  previously  discussed.   Based  on  the  range  of  applicability  of  Plank's 
Equation  of  State  and  comparison  with  the  tables  of  Din,  '  it  is  estimated 
that  temperatures  in  the  vicinity  of  190°K  are  accurate  to  within  .3°K  and 
temperatures  in  the  region  of  l60°K  accurate  to  perhaps  1°  or  2°K.   A  plot 
of  these  points  for  all  three  nozzle  geometries  tested  is  shown  in 
Figure  31. 

Several  characteristics  of  the  condensation  profiles  and  the  local 
conditions  at  point  of  occurrence  may  be  noted. 

First,  the  shape  of  the  profiles  contain  no  sharp  minimum  at  the 

point  of  occurrence  of  condensation  such  as  those  noted  in  the  condensation 

7  23 
of  pure  water  vapor.  '    i.e.  it  is  not  appropriate  to  speak  of  condensa- 
tion shock  in  this  case,  since  the  formation  of  moisture  occurs  along  an 
extended  region. 

Second  it  appears  that  pressure  level  (or  temperature  level,  or  some 
dependent  property)  may  influence  the  shape  of  the  profile,  i.e.  the  rate 
of  formation  of  moisture,  to  some  extent  since  the  shapes  of  the  condensa- 
tion profiles  appear  to  be  more  "spread  out"  for  high  stagnation  pressures 
than  for  low  stagnation  pressures.  See  Figure  27.   These  effects  are  not 
explainable  in  terms  of  gas  dynamics  influences,  and  the  cause  is  not 
apparent . 

Third,  a  very  definite  influence  of  nozzle  geometry  upon  the  point  of 
onset  of  condensation  occurs,  as  seen  in  Figure  31.   The  rapid  rate  of 
nozzle  divergence  and  expansion  of  Nozzle  I  results  in  higher  supersatura- 
tion  ratios  and  relatively  greater  delay  of  condensation.   Progressively 
decreasing  the  rate  of  expansion  through  decrease  of  nozzle  divergence 
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results  in  earlier  condensation.   Comparative  geometry  and  non-condensing 
pressure  profiles  of  the  three  nozzles  may  be  seen  in  Figure  32. 

Fourth,  the  extremely  low  rate  of  divergence  of  Nozzle  II  in  the 
0.80°  case  (subsonic  flow  results  at  about  0.50°  divergence)  results  in 
notable  axial  fluctuation  of  the  pressure  profile  data  (Figure  29)  and 
some  scatter  of  the  point  of  onset  of  occurrence  of  condensation.   (Figure 

31)   It  should  be  remembered  that  differential  changes  in  gas  properties 

2 
are  proportional  to  (l/l-M  )  and  that  at  the  pressure  ratios  present  in 

these  tests  Mach  Number  remains  very  close  to  unity.   (M  -  1.25  at 

P/P   =  .Uo).   Also,  as  heat  of  condensation  is  released,  Mach  Number  is 
o 

driven  even  closer  to  unity.   This  is  coupled  with  the  fact  that  at  low 
angles  of  divergence,  the  value  of  the  negative  pressure  gradient  is 
greatly  decreased,  downstream  of  the  throat,  from  those  values  which 
existed  for  the  other  two  test  nozzles.   This  would  tend  to  create  more 
rapid  boundary  layer  growth  and  greater  departure  from  1-dimensional  flow 
in  the  narrow  angle  case.   It  might  be  expected  that  larger  departures 
from  1-dimensional  flow,  plus  the  magnifying  effects  of  near  unity  Mach 
Numbers  could  contribute  to  larger  pressure  fluctuations. 
2.    Visual  Observations 

Experimental  observations  indicated  that  the  pressure  hump  accompany- 
ing condensation  occurred  at  a  point  coincident  with  the  appearance  of 
the  ice  cloud. 

This  at  first  seemed  at  variance  with  the  expectation  that  pressure 
variation  should  occur  prior  to  the  appearance  of  visible  condensate, 
due  to  the  theoretically  small  sizes  of  forming  nuclei.   (r  -~ 10  cm). 

Yellot's   '    investigation  of  condensing  steam  was  also  accompanied 
by  visual  observation  of  the  water  cloud  coincident  with  the  location  of 
onset  of  condensation  as  determined  by  pressure  measurements.   He 
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reported  that  this  observation  remained  valid  despite  changes  which  resulted 
in  the  axial  shift  of  condensation  in  the  nozzle.  The  same  was  true  in  this 
investigation. 

It  is  also  mentioned  that  the  data  illustrated  a  very  high  degree 
of  consistency  and  reproduceability,  that  illustrated  in  Figures  2k 
through  28  involving  three  separate  sets  of  measurements  over  a  period  of 
six  months,  each  set  being  separated  by  nozzle  dismantlement,  cleaning,  and 
pressure  tap  plate  alterations. 
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IV.   INTERPRETATION  OF  EXPERIMENTAL  RESULTS 

As  discussed  in  Sections  I  and  II,  there  are  reasons  for  supposing  that 
the  process  of  nozzle  condensation  is  one  involving  homogeneous  nucleation 
within  the  super-saturated  vapor.   The  purpose  of  this  section  is  to  discuss 
homogeneous  nucleation  theory  and  the  subsequent  drop  growth  process  and 
apply  these  processes  to  nozzle  expansions  for  the  purpose  of  testing  the 
ability  of  existing  theory  to  interpret  the  experimental  results  of  C0_ 
condensation.  Attention  is  given  to  the  physical  properties  upon  which 
nucleation  and  drop  growth  theory  depend,  and  to  the  uncertainties  which 
exist  in  the  applicable  values  of  these  properties  as  well  as  to  the  theory 
itself. 

A.    Nucleation  Theory 

A  review  of  classical  nucleation  theory,  following  the  development  of 

19 
Frenkel   is  presented  in  Appendix  B.   For  more  comprehensive  review  of 

theory  see  Reference  Numbers  5,6,12,18,38,56,58,65  and  68.  The  discussion 

65 
of  Yang,   accomplished  while  working  under  P.P.  Wegener  at  Yale  University, 

is  particularly  thorough. 

1.   General  Form  of  Nucleation  Rate  Equation 

The  rate  of  production  of  nuclei  in  a  supersaturated  vapor  has  been 

shown  by  various  authors  to  be  of  the  form* 

.     (-AG*/KT)  /_, 

J  =  A»exp  \7) 

where  J  =  number  of  nuclei  produced  per  unit  time  per  unit  volume, 

AG*  =  free  energy  of  formation  of  a  critical  cluster  (a  cluster  in 

unstable  equilibrium  with  its  supersaturated  vapor) 

And  A  is  proportional  to  the  incident  flux  of  vapor  molecules  on  a 

critical-sized  cluster. 

P 


A  ■ Cn  t/is^1  k™*2  (6) 


* 


See  Appendix  B  for  a  review  of  nucleation  theory  development 
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The  value  of  C  varies  somewhat  as  derived  in  the  classical  treatments 
of  Volmer,  '  Becker  and  Doring,   Frenkel,   Barnard,   Yang  ''  and  others. 
Though  variations  of  one  to  two  orders  of  magnitude  occur  over  the  range  of 
values  presented,  the  effect  on  nucleation  rate  is  small,  since  for  only  a 
small  additional  expansion  along  a  nozzle,  the  value  of  J  can  easily  vary  by 
5-10  orders  of  magnitude.   Since  the  change  in  local  conditions  is  small 
for  a  variation  in  nucleation  rate  of  this  sort  and,  in  fact,  is  scarcely 
noticeable  for  1  to  2  order  of  magnitude  variations  in  J,  the  differences 
in  classical  treatments  of  the  theory  are  not  significant. 

In  the  classical  treatments  of  the  theory,  the  nucleus  is  consi- 
dered as  a  stationery  object,  and  the  free  energy  of  formation  is  merely: 

AG  =  W2a   -  gKT  In  |-  ,    . 

00 

where  the  first  term  is  the  work  of  formation  of  surface  area  for  a 
liquid  having  surface  tension  o,  and  the  second  is  the  constant-tempera- 
ture free  energy  change  of  g  molecules  which  change  state  from  vapor  at 
P,T  to  saturated  liquid  at  T,   P   is  the  flat  film  vapor  pressure  of  the 

liquid  at  temperature  T, 

(*) 
The  critical  cluster    is  that  for  which  the  function  AG  is  a  maximum, 

and  results  in 

,   P      2mo 
In  — -  = 


P     pTKTr*  (8) 

°°     L 

Using  the  procedure  presented  by  Frenkel,  for  which  the  value  of  c 
becomes 

C  =  [i^»l1/2  (0) 

the  expression  for  nucleation  rate  is 

T        ,P    ^2   m      ,2c.  1/2        r-UTTor*21  Mnx 

J  =    (KT}      pT  (m}        eXp[      3KT      ]  (10) 

Li 
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As  mentioned,  Yang  has  shown  that  the  results  of  the  authors  previously 
mentioned  yield  expressions  for  classical  nucleation  that  are  much  alike 
(in  some  cases  identical)  and  which  yield  quantitative  predictions  that  are 
not  significantly  different. 

More  recent  authors  have  pointed  out  the  need  for  inclusion  of  addi- 
tional terms  to  the  free  energy  of  formation  of  a  critical  cluster,  which 
account  for  the  fact  that  the  cluster  may  have  translational,  rotational 

and  vibrational  energy  modes „   Kuhrt,  Rodebush  '  Lothe  and  Pound,  '  and 

12  13 
Courtney   '   have  treated  this  problem  but  reach  widely  different  con- 
clusions as  to  the  effect  upon  nucleation  rate.   The  postulated  correction 

factors  for  this  "gassification"  of  the  condensed  clusters  range  from  10 

20 
to  10   ,  depending  on  the  treatment  and  molecular  configuration  chosen. 

i  ft 
More  recently,  Feder,  Russell,  Lothe  and  Pound   have  jointly  made  a 

very  comprehensive  review  of  this  problem,  and  offered  some  new  thoughts  on 
the  subject.  For  water  vapor,  they  conclude  that  the  effects  of  gassifica- 
tion are  equivalent  to  a  correction  factor  of  about  10   ,  but  admit  that 

this  is  not  borne  out  by  cloud  chamber  data.   In  general, such  experimental 

5  23 
data  as  does  exist  tends  to  more  strongly  verify  classical  theory  ' 

than  the  improved  gassified  theory,  despite  strong  arguments  for  the  greater 
theoretical  validity  of  the  latter. 

2.    Evaluation  of  Classical  Nucleation  Rate  Equation 
Essentially  the  same  physical  properties  appear  both  in  the  exponen- 
tial and  pre-exponential  terms  of  the  nucleation  rate  equation.   Of  course, 
uncertainties  in  these  properties  will  effect  the  calculated  value  of 
nucleation  rate  chiefly  through  the  exponential.   By  substitution  for  the 
value  of  r*  it  is  seen  that  the  exponential  is 

r-UTTar*2 1  r-UiT    ,        2m        \  2  / a    \  3 1 

expL      3KT      J   =  expL    3      ^p    InP/pJ    ^KT;    J  (ll) 
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In  evaluating  this  term  correctly  it  is  clear  that  four  physical 
properties  of  the  vapor  must  be  known:  l)the  liquid  drop  surface  tension 
2) local  temperature  3) liquid  drop  density  and  h) local  supersaturation 
ratio. 

Much  attention  has  been  given  to  the  applicable  values  of  surface 
energy  for  small  droplets.   However,  here  again  there  is  little  agreement 
as  to  what  corrections  should  be  made  to  flat  film  values  to  take  account 
of  the  effect  of  droplet  curvature,  and  the  sign  of  the  correction  is  itself 
a  subject  of  controversy. 

Effects  of  curvature  are  considered  by  some  to  cause  a  surface  tension 
decreasing  with  drop  radius  such  as  in  the  Kirkwood-Buff  Equation. 

0  =  (1  +  26/r)  (12) 

where  6  is  a  length,  nearly  constant,  lying  between  .25  and  0.6  of  the 

molecular  radius. 

56 
In  addition,  Stever  and  Rathbun   have  calculated  a  correction  to  take 

account  of  the  radius  dependency  of  surface  tension  which  also  results  in  a 

decrease  for  small  radii.   Compared  with  uncorrected  surface  tension  these 

corrections  lead  to  lower  energies  of  formation  of  nuclei  and  hence  higher 

nucleation  rates  for  the  same  supersaturated  conditions. 

kk 
However,  Oriani  and  Sundquist,   predict  an  increase  in  surface  energy 

of  about  25%  in  the  case  of  water,  using  an  analysis  based  on  the  breaking 

of  intermolecular  bonds.   A  bonding  model  used  by  Benson  and  Shuttleworth, 

however,  predicts  a  15%  reduction  in  surface  energy  of  correspondingly  small 

o 

water  ice  crystals.   Bogdonoff  and  Lees,   also  using  considerations  of  bind- 
ing energies  of  sodium  and  of  nitrogen  obtained  a  correction  that  results 
in  larger  critical  clusters  and  also  larger  energies  of  formation,  hence 
lower  nucleation  rates  from  the  uncorrected  case.   For  lack  of  a  more 
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definite  indication  of  the  effects  of  small  droplet  sizes,  most  investiga- 
tors use  the  flat  film  value  of  surface  energy,  and  estimate  the  effects 
on  the  theoretical  predictions  of  departure  from  it. 

An  additional  uncertainty  in  surface  energy  occurs  with  CO  due  to 

the  fact  that,  at  pressures  below  its  ice  point,  75.1  psia,  it  apparently 

9 
condenses  in  the  form  of  ice  crystals  rather  than  liquid  droplets.   Bondi 

has  shown  that  for  metal  crystals,  which  involve  valence  bonding,  increases 

in  surface  energy  from  the  liquid  value  are  closely  proportional  to  the 

heat  of  fusion  and,  depending  also  on  the  crystal  planes  involved,  amount 

to  increases  of  about  25%  for  those  metals  which  he  discusses.   He  further 

presumes  that  for  molecular  solids  one  might  expect  increases  of  from 

approximately  1/8  to  1/3  these  amounts. 

In  the  case  of  some  gases  there  is  doubt  as  to  the  proper  value  of 

droplet  density  to  be  applied.   Though  density  does  not  vary  too  rapidly 

with  temperature,  there  is  doubt,  in  the  case  of  CO  ,  as  to  whether  solid 

or  liquid  density  should  be  used.   Although  the  CO  may  condense  as  ice 

particles,  there  is  some  argument  for  use  of  liquid  density  in  order  to  be 

consistent  with  the  spherical  geometry  and  liquid  surface  energy  applied  to 

the  liquid  drop  model,  and  in  the  absence  of  any  knowledge  of  CO   crystal 

surface  energies.   Also,  studies  of  Dorsch  and  Hacker   with  very  small 

water  droplets  indicated  that  the  spontaneous  freezing  temperature 

decreases  rapidly  with  decreasing  droplet  size.  (Figure  33)  This  might  tend 

to  justify  doubt  on  the  reality  of  attempting  to  distinguish  the  difference 

between  liquid  and  solid  states  at  droplet  sizes  of  critical  clusters.   [For 

p  7 

CO  ,  typical  critical  radii  are  of  7.5  x  10~  cm  to  10~  cm  (about  30  -  50 

molecules)  at  supersaturation  ratios  of  the  order  of  10,  where  noticeable 

17 
quantities  of  moisture  are  condensed.   Typical  nucleation  rates  are  10 

21  3 

to  10   nuclei/cm  -  sec] 
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Figure  3^  illustrates  existing  data  on  vapor  pressures  of  CO   in  the  region 
of  experimental  condensation,  and  the  computational  approximations  used  in 
the  condensation  calculations.  Accuracy  of  the  vapor  pressure  data  and  the 
approximation  below  temperatures  of  about  l80°K  is  uncertain,  but  is  probably 
within  5%   for  higher  temperatures.   Supersaturation  ratio  may  be  determined 
within  these  limits  of  accuracy  for  a  given  local  pressure  if  the  tempera- 
ture is  also  known. 

Surface  tension  is  considered  to  be  described  for  flat  film  liquid  by 
the  function  o^  =  5»00  -  .193T(°C)  dynes. /cm,  based  on  data  of  References  21, 

27,  and  50,  for  which  a  linear  dependence  of  o«>   on  T  down  to  220°K  is 

T   125 
justified.   An  alternative  function  o  =  75(l-m~")  "'  (Reference  27)  yields 

c 

values  of  aoo  which  vary  by  only  2%   to  9%9   in  the  200°K  to  l60°K  region  of 

interest,  from  those  of  the  linear  function  above. 

Knowledge  of  local  temperature,  previously  discussed,  is  important  in 
evaluating  the  nucleation  rate  equation,  not  only  because  of  its  presence  as  a 
cul  i.c  iii  bhe   .;;  no  I  ,  but  because  bhe  other  variables  are  also  dependent  in 
value  upon  temperature.   Effects  of  uncertainty  in  temperature  do  not  cancel 
appreciably  in  their  influence  upon  theoretical  nucleation  rate.   At  200°K, 
a  1%   error  in  temperature  will  result  in  about  a  2%   error  of  the  opposite 
sign  in  surface  tension  and  hence  about  a  3%   error  in  the  ratio  o/T  for 
CO  .  These  effects  were  illustrated  in  the  earlier  discussion  of  the 
potential  uncertainty  of  local  temperature  unless  an  equation  of  state  is 
used  to  determine  the  isentrope  for  CO  .   It  is  possible  that  some  past 
applications  of  theory  to  experimental  condensation  data  have  erred  through 
failure  to  consider  these  real  gas  effects. 

Departures  from  perfect  gas  behavior  affect  the  evaluation  of  the  free 
energy  required  for  the  formation  of  a  critical  cluster,  and  also  the  drop 
growth  relations  dealt  with  in  the  following  section,  since  both  nucleation 
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and  drop  growth  developments  are  based  on  perfect  gas  assumptions.,   This 
problem  is  dealt  with  in  Section  IVD,  following  the  discussion  of  drop 
growth  theory,  and  the  application  of  nucleation  and  drop  growth  to  condensa- 
tion in  nozzles. 

B.   Drop  Growth  Theory 

Subsequent  to  the  production  of  nuclei,  additional  formation  of 
moisture  will  occur  through  the  growth  of  the  clusters  due  to  vapor  impinge- 
ment.  Nuclei  may  by  this  process  grow  to  radii  of  up  to  about  an  order  of 
magnitude  greater  than  their  original  size,  depending  upon  two  important 
drop  growth  parameters  introduced  below.   In  the  theoretical  treatment,  the 
distinction  between  nucleation  and  commencement  of  drop  growth  treatment  is 
somewhat  artificial,  as  it  is  necessary  to  establish  some  size  of  cluster 
at  which  the  nucleated  droplet  exists  and  may  commence  growing.   Becker 
and  Doring  have  shown  that  at  a  cluster  radius  1.3  times  critical,  the 
probability  of  clusters  continuing  to  grow  reaches  better  than  99%,  and 
hence  it  is  considered  that  this  is  a  reasonable  starting  point  for  the 
application  of  drop  growth  equations.   (Note:   Computations  showed  that 
results  were  relatively  insensitive  to  this  starting  radius.) 

The  growth  of  a  droplet  depends  predominantly  on  two  factors:  l)The 
rate  of  condensing  mass  flux  2)The  rate  of  droplet  energy  transfer  to  the 
vapor  environment. 

It  would  appear  that  the  first  factor  is  itself  a  description  of  the 
drop  growth  process,  but  in  actuality,  it  is  greatly  dependent  upon  not 
only  local  vapor  conditions,  but  also  the  difference  between  droplet  and 
vapor  conditions.   Both  of  these  are  strongly  influenced  by  the  second 
factor. 
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Two  coefficients  are  defined  to  describe  this  behavior: 

£  =  Condensation  Coefficient.   The  fraction  of  impinging  vapor 

molecules  which  are  absorbed  at  the  drop  surface, 
a  =  Thermal  Accommodation  Coefficient.   The  fractional  temperature 

change  which  takes  place  in  the  reflected  portion  of  incident 

molecules. 

«   =  m m  (13) 

D~ 

Temperature  variation  across  the  droplet  is  assumed  to  be  negligible, 
and  the  droplet  is  assumed  to  be  spherical. 

Based  on  mass  and  energy  balance  considerations  presented  and  developed 
in  Appendix  C  the  following  drop  growth  equations  are  obtained. 


Conservation  of  Mass 

pLUnr  dF  =  U*r  (T2^r7l/2  "  T2^yi/2} 
Conservation  of  Energy 

U/3lTr3pLc£dT  +  kn*7£   PLUfD  =  6(2RT)  "  ^D(2RTD} 
-  (1  -C)62R[T  +  a(TD  -  T)]  (15) 

These  equations  may  be  reduced,  respectively,  to 

Li 

o  =  2C(i  -jpffc)  -  2(1  -oot  (^- 1)  -  e(i  -f11)  (7-^-1)  <x  - 1)    (it) 

where  A   =  I-=—    M& 

a  _Z ft        _  PD  . 

P   "    (2^RT)l/2  %  ~    (2*RTDH/2 

And  P   is  the  equilibrium  vapor  pressure  for  a  spherical  liquid  drop  of 

temperature  T^  in  an  environment  of  temperature  T^. 

D  on  D 

(      2°      ) 

PD  =  P^  (TD)  expVPLRTDr  (18) 

In  the  further  application  of  these  equations  to  streamwise  growth  of 
drops  in  a  nozzle,  additional  assumptions  are  made.   First,  for  an 
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acceptably  simple  computation  scheme  it  is  convenient  to  treat  droplets  of 
uniform  size,  where  the  radius  is  an  average  based  either  on  surface  area 
or  volume.   Ignoring  the  effects  of  coagulation,  there  is  evidence  that  the 
actual  distribution  of  drop  sizes  is  narrow,  and  that  the  average  droplet 
treatment  does  not  introduce  appreciable  error. 

One  possible  error  that  may  occur  as  a  result  of  the  above  assumption 
becomes  likely  when  the  supersaturation  ratio  P/P^,,  is  reduced  so  rapidly 
as  a  result  of  condensation,  that  droplets,  which  have  recently  nucleated 
and  begun  to  grow,  will  suddenly  be  of  size  less  than  the  critical  cluster 
size.   In  this  case,  considerable  re-evaporation  would  occur  among  the 
smaller  droplets  of  the  actual  size  distribution,  and  the  treatment  of  an 
"average"  droplet  does  not  account  for  this. 

Relative  motion  between  droplets  and  the  gas  stream  is  ignored,  and 
though  this  is  generally  accepted,  it  may  be  possible  that  the  same 
translational  modes  of  energy  required  by  recent  authors  in  nucleation 
theory  might  also  affect  to  some  degree  such  things  as  incident  mass  and 
energy  fluxes  upon  newly  formed  nuclei. 

Values  of  a  and  £  have  been  presented  by  numerous  investigators  for 
several  vapors,  excluding  CO  .   Although  earlier  presented  values  of  £  were 
fairly  low,  there  has  been  more  recent  substantiation  of  values  close  to 
1.00  for  most  materials,  including  liquid  metals,  the  previous  errors 
having  been  caused  by  neglect  of  surface  temperature  depressions  and  the 
effects  upon  evaporation  flux.   However,  there  is  also  indication  that 

molecules  of  higher  dipole  moment  tend  to  have  lower  values  of  £,  as  with 

12  UQ 

liquid  water  for  which  Alty  and  MacKay,  '  and  Prueger   present  values  of 

less  than  £  =  o0U„   The  influence  of  crystallization  upon  K  may  also  be 

large,  as  Tschudin  presents  a  value  close  to  unity  for  water  ice.   Effects 
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of  surface  contamination  and  surface  diffusion  barriers  (to  recently- 
evaporated  particles  due  to  insufficient  vacuum  and  mean  free  path  lengths,) 
are  also  thought  to  play  a  role  in  experimentation  which  tends  to  lower 
apparent  values  of  C.  No  experimental  values  of  £  have  been  measured  for 
CO  ,  which  has  a  dipole  moment  of  about  1.5  compared  with  80.0  for  that  of 
water. 

The  thermal  accommodation  coefficient  for  vapors  and  their  liquids 
is  generally  accepted  to  be  around  unity.   However,  values  may  be  much 
lower  between  vapors  and  solids,  as  Wachman   indicates  for  gases  on  clean 
metal  surfaces.   For  CO  ,  the  value  of  a  for  the  vapor  and  its  solid 
crystal  is  unknown,  and  there  is  uncertainty  as  to  which  of  the  two  situa- 
tions above  is  the  most  applicable  description. 

Departures  from  perfect  gas  behavior  also  influence  drop  growth. 
Again,  correct  values  of  local  temperature  and  supersaturation  ratio  will 
be  influenced,  and  the  perfect  gas  assumptions  used  in  treating  mass  and 
energy  balances  leading  to  equations  (l6)  and  (IT)  may  require  modification. 

In  this  study,  considerable  attention  has  been  given  to  the  behavior 
and  influence  of  the  above  mentioned  factors. 

C„   Application  of  Nucleation  and  Drop  Growth  Equations 

The  predicted  effects  of  moisture  formation  during  a  nozzle  expansion 
are  determined  by  applying  the  time  rate  release  of  heat  to  the  gas  dynamics 
of  the  related  expansion.   By  considering  the  process  with  respect  to 
variable  nozzle  position,  rather  than  time,  a  result  is  obtained  which  can 
be  directly  compared  with  experiment. 

Interferometric  measurements  of  density  variations  in  the  nozzle 
verified  that  the  flow  between  boundary  layers  is  in  fact  one-dimensional 
and  allows  the  use  of  the  one-dimensional  gas  dynamics  relations; 
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Continuity 

pAu  / , „ » 

Momentum 

-Adp  =  mdu  (20) 

Energy 

C  dT  +  udu  =  h   dy  (21) 

P  fg 

where  m  is  the  mass  rate  of  flow  through  the  nozzle,  V  the  mass  fraction 

of  condensed  moisture,  u  is  the  streamwise  velocity  of  the  vapor,  and  p, 

P,  and  T  are  the  vapor  density,  pressure  and  temperature,  respectively. 

A  is  the  effective  area  of  flow  between  boundary  layer  thickness,  and  is 

determined  experimentally  for  non-condensing  flow. 

The  rate  of  moisture  formation  is  determined  as  follows  and  assumes 

the  droplet  velocities  to  be  those  of  the  vapor.   A  nucleus  of  radius  r 

formed  at  some  point  x  will  grow  to  a  larger  size  at  some  future  time  and 


distance  x, 


+  r   I  dr 


dx  (22) 

u    x  u  dt 
o 

and  the  surface  area  will  be 

a  =  Uirfr  +  '   -~  dx]2  (23) 

o    x  u  dt 
o 

dr 
For  a  rate  of  growth  —  predicted  by  drop  growth,  the  increase  in 

liquid  mass  per  droplet  will  be 

-~dx]2|£±dx  (2U) 

L  L  o  x     u  dt  dt   u 

o 

And  the   fractional  moisture  change  that  occurs   at  x  for  the 

dx 
dn     =  J(x    )A(x    )[ — —   ]dx     particles   formed  in  the  volume  A(x    )dx      is 
o  oou  o  -,  o,onJ 

°  £^2.  r         fX     ^—  dxl2—  — 

^pPLadr  =  J(x0)A(x0)   uq     dx0PLUirLro'>'Jxo  u  dt        J    dt  u 

%  "  (P  A  dx  )/(l-M  )  (p  Adx)/(1  -  u  ) 

o  o  o      o  ^o  o  o        o 

(25) 
Inclusion  of  the  drop  growth  on  all  previously  formed  nuclei  necessi- 
tates integrating  the  right  hand  term  over  all  x  .   Noting  that 


32. 


m  =  p  A  u  /l  -  u   and  considering  also  the  moisture  produced  through 
o  o  o      o 

additional  nucleation  at  x,  there  results 

111  =  f  J(x  )(A(xa)dxa)pLU]L  [r  +!Xl^Ldx]2^dx 
dx     o    o     u       m    o    x  u  dt      dt   o 

+  l1--*3-  i,(i)A(x) 
3     m 

WL   rf    /  r      1  dr        \?  T/       \./       \    1  dr   , 

=  -t-^  [J    (r     +  J      -  —  dx )/  J  x     A(x       -  —  lx 
m          ^  o     o        '  x  u  dt  o  oudto 

o 

+   1/3  r3  J(x)A(x)]  (26) 

For  nuclei  generally  only  slightly  larger  than  critical  size,  the 

23 

theoretical  drop  growth  rate  is  essentially  independent  of  size.    Hence 

dr 

—  is  assumed  to  have  a  value  equal  to  that  for  large  drops.   Also,  since 

it  is  practically  impossible  to  keep  account  of  the  radii  of  each  group 

of  growing  nuclei  which  formed  at  each  section  dx,  it  is  assumed  that  the 

droplets  can  be  represented  by  the  "surface-area  averaged"  radius  of 

droplets.   As  previously  mentioned,  volume  averaging  has  little  influence 

on  the  results.   In  addition,  Lundgren  has  given  experimental  evidence 

that  the  actual  size  distribution  of  particles  condensed  in  a  nozzle 

expansion  is  quite  narrow. 

The  foregoing  equations  are  reduced  to  the  following  form  for  use 

with  a  Runge-Kutta-Merson  integration  procedure  along  the  expansion  path. 

(The  detailed  listing  of  computation  procedure  may  be  seen  in  Appendix  H). 

dYj    =   8ttJA 

dY2   =   Yl   f  +  r*dYl 
dY,    =   Y?   |f  +  JA(U7rr*2) 

«5  =  I  -  p[<a  -  pM*  -  ££i  '  [1  - (1  -  ^)(Y^ +  w)] 

dYr>    =  £     =  T[Ad,«    +  Y~-A(l   _„)   |P] 
'        dx  l"  M2  P  ?TJ 
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where 

rg  P 

y  =   c  /c 

P  v 
M  =  u//yRT 

Initial  values  of  the  stagnation  conditions,  physical  properties  of  the 
vapor  and  its  liquid,  values  of  a  and  £ ,  and  calculation  starting  point  are 
given  for  each  computation .   Effects  of  contamination  may  also  be  considered 
by  giving  non-zero  values  to  Y(l),  Y(2)  and  Y(3),  based  on  an  initial  concen- 
tration and  size  of  contaminant.   Effects  of  possible  droplet  energy  modes 
upon  free  energy  of  cluster  formation  ("gassified"  behavior  of  droplets  - 
page  23)  may  be  considered  by  multiplying  nucleation  rate  by  constant  arbit- 
rary factors,  equivalent  to  the  values  proposed  by  various  authors. 

The  various  sub-routines  of  the  listed  computation  scheme  accomplish 
the  determination  of  such  things  as  the  change  in  nozzle  area,  vapor  and 
liquid  properties,  super saturation  ratios  and  solution  to  the  drop  growth 
equations  when  step-wise  changes  are  made  along  the  stream  of  the  expansion. 
An  error  estimation  is  made  for  each  computation  and  the  size  of  the  step, 
dx,  is  adjusted  to  reduce  the  error  bounds  to  within  the  desired  limits. 

D.    Real  Gas  Effects 

The  preceeding  developments  of  nucleation  and  drop  growth  theory  and 
application  to  nozzle  condensation  computations  assume  perfect  gas  behavior 

during  the  expansion.   However,  as  earlier  discussion  illustrated,  experim- 

Pv 
ental  measurements  were  being  made  for  which  values  of  —  which  were  consid- 

KI 

erably  less  than  unity,,   Consequently,  the  computational  scheme  was  modified 
so  that  the  effects  of  departure  from  perfect  gas  behavior  could  also  be 
considered. 

A  feature  of  an  isentropic  expansion  which  facilitates  fairly  simple 
corrections  to  the  gas  dynamics,  drop  growth  and  nucleation  equations  is  the 
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fact  that  for  this  particular  thermodynamic  process,  the  compressibility- 
factor  remains  fairly  constant,  changes  of  less  than  2%   in  the  value  of  z 
generally  occurring  for  expansion  pressure  ratios  of  10  or  more.   This  was 
found  to  be  true  for  several  other  gases,  including  ammonia,  and  for  CO 
is  valid  throughout  the  superheated  and  supersaturated  regions  for  which 
compressibility  factors  are  in  excess  of  about  .88.   Using  this  fact,  it 
may  be  shown  (Appendix  D)  that  the  following  relations  are  approximately 
true  for  an  isentropic  expansion. 

Pv  ,   \ 

—  =  z  ^constant  (27) 

*'-  i ;  (,y-d(i-z)  -  f  <?»  <28> 

R'  =  p-=     zR  (29) 

m 

where  y1  and  R'  are  the  effective  ratio  of  specific  heats  and  gas 
constant,  respectively,  that  are  applicable.   It  should  be  noted  here,  that 
the  nucleation  rate  is  effected  in  the  following  manner; 


J  <   exp  (         ) 


i    „?  2mo 


* 


3KT    >'  pTzKT  £n  P/P 

L  °° 

due  to  the  change  in  free  energy  required  to  form  a  critical  cluster. 

To  apply  these  effects  to  the  perfect  gas  relations  of  the  condensation 
computations,  two  additional  steps  were  necessary.   First,  modifications  which 
allowed  for  use  of  a  variable  ratio  of  specific  heats  (Y  =  Y(T))  were  made. 
Second,  each  set  of  stagnation  conditions  which  was  to  be  placed  in  the 
condensation  scheme  first  had  its  "real"  isentrope,  z,  and  effective  isentro- 
pic exponent  Y'(.T)  determined,  and  from  these  the  Y(T)  which  corresponded 
to  z  and  Y1,   Computation  input  was  Y(T)  and  correction  (compressibility 
factor)  z.   This  resulted  in  an  isentrope,  providing  the  correct  values  of 
P  vs.  T  being  employed  in  each  condensation  calculation,  plus  the  use  of 
effective  values  of  y   and  R  in  the  drop  growth  and  nucleation  calculations. 
Additional  details  of  this  procedure  are  presented  in  Appendix  D. 
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E.    Results  of  Application  of  Theory  to  Condensation  of  CO? 
1.    Comparison  of  Theory  and  Experimental  Data 

Data  shown  in  Table  lU  and  Fig.  31  illustrate  the  local  conditions  at 
the  occurence  of  condensation  for  each  of  the  experimental  tests  made. 
Local  temperature  was  determined  through  computation  of  the  isentrope  for 
each  test  case,  using  the  procedure  discussed  on  page  13  based  on  Plank's 
Equation  of  State  for  C02 .   Reference  to  Figure  35  illustrates  that  .  1% 
presence  of  moisture  is  a  convenient  and  reasonable  criteria  for  theoretical 
onset  of  condensation.   Comparing  the  data  and  the  predictions  of  theory- 
involved  two  main  comparisons;  l)  The  location  of  the  onset  of  condensation, 

defined  as  the  theoretical  occurrence  of  .1%   moisture,  and  2)  The  shape  of 
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the  pressure  profiles.   For  lack  of  better  knowledge,  <*,  K   and  —  were 
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originally  assumed  to  be  unity.   Specific  comparisons  were  made  of  the  pred- 
ictions of  l)  Classical  nacleation  theory  2)  Droplet  "gassification" 
corrections  to  the  theory   3)  "Real"  gas  vs.  psrfect  gas  treatment  of  the 
classical  theory  (to  determine  whether  or  not  a  simplified  approach,  contain- 
ing known  errors,  might  yet  yield  a  better  prediction  of  observed  behavior.) 
U)  Use  of  liquid  vs.  solid  droplet  density  and  5)  the  predictions  of  the 
theory  in  nozzles  of  different  geometry.   Figures  36  through  39  illustrate 
these  comparisons  and  support  the  following  conclusions? 

1)  All  predictions  based  on  values  of  c/a^,  a  and  £  of  unity 
predict  condensation  at  considerably  lower  supersaturation  ratios  than  those 
observed  experimentally. 

2)  Real  Gas  treatment  of  the  classical  theory  yields  a  closer 
prediction  than  perfect  gas  treatment,  and  the  "gassification"  corrections 
are  in  the  wrong  direction.   (Fig.  36) 

3)  Use  of  liquid  rather  than  solid  droplet  density  yields  a 
prediction  nearer  to  experimental  results,  lacking  any  better  knowledge  of 
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what  corrections  might  be  applied  to  surface  energies  of  small  C02  crystals. 
(Fig,  37) 

M   The  unmodified  theory  does  extremely  poorly  in  predicting  any 
effects  of  nozzle  geometry.   (Fig.  38) 

5)   The  shape  of  condensation  profiles  predicted  by  values  of  « 
and  £  of  unity  are  in  extreme  disagreement  with  those  observed  experimentally. 
(Figs.  U5  and  U6) 

2.   Influence  of  Nucleation  and  Drop  Growth  Theory  on 
Theoretical  Predictions. 
It   is  concluded  that  homogeneous  nucleation  and  drop  growth  theory  fail 
in  their  ability  to  predict  observed  condensation  behavior  of  CO  ,  unless 
some  modification  in  values  of  the  nucleation  and  drop  growth  parameters  a, 
p  ,  a,  and  £  can  be  independently  established.   Toward  better  understanding 
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of  how  these  parameters  influence  predicted  condensation  phenomena,  and  with 
the  idea  in  mind  of  finding  corrections  for  the  prediction  of  CO  condensa- 
tion, a  series  of  calculations  based  on  arbitrary  variations  of  o,  pT,  a  and 

Li 

£  were  made. 

Figures  39  and  UO  show  typical  pressure  profiles  for  condensation  of 
C0_,  for  various  condensation  and  thermal  accommodation  coefficients.   Shown 
in  the  lower  portion  of  the  figures  is  the  nucleation  rate,  indicating  how 
nucleation  rate  remains  high  for  longer  duration  when  drop  growth  is  retarded 
This  results  in  condensation  clouds  of  greater  particle  density  and  smaller 
particle  size  than  in  the  case  of  ready  drop  growth.  In  addition,  smaller 
total  deflection  in  the  static  pressures  (from  non-condensing),  indicating 
lesser  amounts  of  moisture,  exist  for  the  downstream  portions  of  profiles 
having  lower  drop  growth  parameters. 

An  important  observation  in  these  figures  is  that,  beyond  a  certain 
reduction  in  the  drop  growth  parameters,  appreciable  proportions  of  moisture 
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may  be  present  primarily  from  nucleation.   This  is  manifested  in  the  form 
of  a  so-called  "ramp"  deflection  in  the  pressure  profile  to  which  additional 
drop  growth  effects  are  super-imposed.   Thus,  for  example  in  Figure  39,  for 
an  assumed  a  of  .10,  variations  in  £  from  1.0  to  .3  and  then  to  .1  have 
progressively  less  effect  on  the  point  of  departure  of  the  condensation  curve, 
and  for  lower  values  of  £,  the  nucleation  "ramp"  becomes  increasingly  appar- 
ent.  The  same  behavior  occurs  for  decreasing  values  of  a,  (shown  for  an 
assumed  value  of  £  =  .  0U)  in  Figure  Uo.   It  would  be  possible  to  obtain  a 
deflection  curve  involving  only  nucleation,  if  drop  growth  were  set  at  zero. 

Figure  4l  illustrates  the  influence  of  increasing  or  decreasing  the 
nucleation  rate  by  changing  the  free  energy  of  formation  of  nuclei.   Indicated 
are  variations  in  the  surface  energy  from  flat  film  liquid  value,  though 
similar  changes  may  be  introduced  by  varying  the  droplet  density  p  .  The 
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important  parameter  in  the  nucleation  rate  exponent  is  (p3/pT  ),  since,  with 
knowledge  of  temperature  and  supersaturation  ratio,  the  uncertainties  of 
surface  energy  corrections,  small  cluster  densities,  and  non-spherical  geom- 
etry may  all  be  contained  in  this  term. 

Two  interesting  features  may  be  noted  in  Figure  kl.      First  is  the  effects 
of  higher  downstream  velocities,  which  tend  to  "wash  out"  the  curves  further 
downstream,  and  second,  the  important  observation  that,  for  a  and  £  of  unity, 
condensation  results  in  a  return  to  near-equilibrium  conditions,  seen  by  the 
fact  that  condensation  commencing  at  any  point  in  the  nozzle  results  in  the 
same  pressure  deflection  immediately  further  downstream, 
a)   Onset  of  Condensation 

As  the  previous  figures  suggest,  influence  of  the  drop  growth  coeffic- 
ients upon  the  point  of  onset  of  condensation  is  relatively  small  compared 
to  that  caused  by  changes  in  surface  tension.   However,  uncertainties  in  a 
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and  £  still  have  appreciable  effect  upon  values  of  surface  energy  necessary 
to  bring  the  predictions  of  theory  into  agreement  with  experiment .   Figures 
1|2,  ^3,  and  JtU  show  the  corrections  to  flat  film  surface  tension  necessary 
to  provide  agreement  between  theory  and  experiment  (based  on  fitting  the 
data  of  nozzle  11(2.20°)).   1*  the  uncertainty  in  the  values  of  a  and  £  for 
CO  crystals,  is  given  limits  of  from  perhaps  .01  to  1.0  for  each,  then  the 
respective  surface  energy  corrections  necessary  may  range  from  about  2Qf% 
to  38%  increases  in  the  flat  film  surface  tension.   For  the  absolute  lower 
bound  of  uncertainty  on  a  and  £  of  zero,  the  necessary  correction  in  surface 
energy  is  reduced  to  +17%  . 

~L    is  noted  in  these  figures  that  the  theoretical  effects  of  nozzle 
geometry  are  increased  to  some  extent  by  use  of  small,  but  non-zero  values 
of  the  drop  growth  parameters,  but  that  the  theory  in  all  cases  predicts  less 
effect  than  that  recorded  experimentally. 
b )  Profile  Shapes 

Figure  U5  illustrates  the  pressure  profiles  for  a  single  expansion 
corresponding  to  the  same  values  of  a  and  £,  and  corresponding  o  determined 
above,  and  indicates  an  additional  distinct  failing  of  the  theory  in  its 
description  of  the  early  stages  of  condensation.   1  is  apparently  impossible 
to  reconcile  the  experimental  downstream  pressure  deflections,  which  indicate 
a  considerable  degree  of  super saturation,  or  departure  from  equilibrium,  with 
the  apparent  rates  of  formation  of  moisture  immediately  after  onset.   When 
drop  growth  parameters  are  adjusted  to  bring  theoretical  profiles  into 
agreement  regarding  the  amount  of  downstream  pressure  deflection,  then  the 
upstream  theoretical  rate  of  formation  of  moisture  lags  considerably  behind 
the  experimental  results. 

This  problem  may  result  from  some  of  the  approximations  discussed  in 
the  numerical  calculations,  or  may  in  fact  indicate  some  physical  "arrest" 
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of  the  condensation  process,  perhaps  due  to  decreasing  mass  or  thermal 
accommodation  characteristics  of  crystals  which  grow  rapidly  colder  during 
their  passage  downstream. 

c)   Downstream  Moisture  Formation 
For  purposes  of  calculating  the  amounts  of  carbon  dioxide  that  will 
form,  or  the  downstream  pressure  variations  which  will  occur  during  a  nozzle 
expansion  of  CO  ,  the  following  parameters  were  determined,  having  arbitrar- 
ily set  thermal  accommodation  coefficient,  a,  at  unity.   There  results,  as 
a  best  description  of  the  three  nozzles,  the  profiles  shown  in  Figures  ^6, 

1*7,  hQ   and  Ua  for  a  =  1.0,  £  =  .00008,  pT  =  1.18  gm/cm3 ,  p/p  =  1.18.   These 
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drop  growth  parameters  have  no  physical  significance,  and  are  presented  only 
as  an  aid  in  determining  approximate  amounts  of  CO  moisture  formation  sub- 
sequent to  the  occurrence  of  condensation  in  a  nozzle  expansion. 

The  best  compromise  to  profile  shape  soon  after  onset  has  been  a  recipe 
which  yields  theoretical  profiles  lagging  behind  and  below  the  experimental 
data  early  in  the  condensation,  and  which  then  pass  slightly  to  the  other 
side  of  the  data  by  about  the  same  amount  for  the  remainder  of  the  profile. 

The  condensation  data  of  Nozzle  II  (0.80)  indicates  the  presence  of  >«r 
some  cyclic  axial  variation  in  the  pressure  ratios  of  all  profiles.  This 
suggests  that  at  this  narrow  angle,  some  departure  from  2-dimensional  flow 
may  be  interfering  with  the  condensation  process.   Friction,  or  compression 
shocks  may  be  creating  these  fluctuations,  and  would  also  result  in  decreased 
supers aturat ion  ratios  in  this  region,  offering  a  possible  explanation  as  to 
why  considerably  less  moisture  is  forming  than  the  recipe  predicts  for  that 
nozzle. 
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3 •   Influence  of  Contamination 

Computations  attempting  to  determine  the  influence  of  contamination  of 

vapor  in  nozzles  were  carried  out  at  intervals  during  this  investigation. 

Original  calculations  were  based  on  perfect  gas  assumptions  and  values  of  a, 

£,  and  -  of  unity,  in  an  attempt  to  establish  whether  or  not  condensation 
c 

CO 

could  possibly  be  due  to  reasonable  amounts  of  contamination  present  in  the 
test  vapor.   Such  contaminants  might  be  present  in  the  form  of  dust,  charged 
particles,  or  the  result  of  an  earlier  condensation  of  some  other  vapor. 
Contamina.tion  of  this  sort  would  be  expected  to  cause  earlier  condensation 

than  would  occur  if  the  test  vapor  were  perfectly  homogeneous. 

27 

Junge  presents  data  for  typical  size  density  distributions  of  contam- 
ination in  atmosphere.   Upper  limits  of  particle  concentration  in  dirty 
atmosphere  of  a  large  city  reach  values  of  5  x  105  cm"'   for  particles  of 
size  less  than  10  5  cm  (.1  micron),  about  10?  cm   for  sizes  from  .1  to  1.0 
microns  and  1.0  cm    for  particles  exceeding  one  micron  in  size. 

Typical  profiles  of  supersaturation  obtained  for  increasing  concen- 
trations of  contamination  with  microscopic  particles  are  shown  in  Fig.  50. 
Beyond  a  limiting  concentration  of  particles,  (depending  on  size)  the  super- 
saturation  ratio  obtains  lesser  and  lesser  values,  until  the  flow  reaches 
a  state  where  it  is  close  to  equilibrium  throughout  the  expansion,  due  to 
the  r<=ady  formation  of  moisture  on  contaminant  particles. 

Subsequent  to  the  establishment  of  the  "recipe"  described  in  the 
previous  section,  additional  computations  for  both  microscopic  and  macro- 
scopic contamination  were  made. 

The  results  are  shown  graphically  in  Figures  51  and  52,  and  are  briefly 
summarized: 

l)  Micron-sized  particles  in  excess  of  108/  cm3  are  required  to 
influence  the  theoretical  condensation  of  C02.  This  is  at  least  5  orders  of 
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magnitude  larger  than  upper  atmospheric  concentrations  of  particles  of  this 
size. 

2)  Microscopic  particles,  such  as  might  occur  from  the  preconden- 

_7 
sation  of  another  vapor,  of  size  the  order  of  10  cm  are  required  in  excess 

of  1012/cm  to  influence  theoretical  condensation  of  C02 . 

3)  No  appreciable  difference  between  empirically  fitted  real  and 
perfect  gas  theories  results. 

h)     The  difference  between  amounts  of  contamination  necessary  to 
influence  condensation  in  the  rapid  expansion  of  Nozzle  I  and  the  slow 
expansion  of  Nozzle  II  (0.80  )  is  approximately  one  order  of  magnitude, 
independent  of  particle  size. 

5)  The  empirically  fitted  recipe  predicts  necessary  concentrations 
of  particles  roughly  two  orders  of  magnitude  larger  than  that  predicted  using 
values  of  a,  £,  and  o/0^  of  unity. 
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When  the  curves  of  Fig  52  are  replotted  in  terms  of 
two  curves  tend  to  move  together,  sinre  the  first  curve  to  drop  downward 
represents  larger  volume  contamination.  This  results  in  fairly  localized 
values  of  supersaturation  ratio  at  occurrence  of  condensation  as  a  function 
of  volume   contamination,  regardless  of  particle  size. 

Applying  the  volume  amounts  of  contamination  of  CO   in  nitrogen  presen- 
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ted  by  Willmarth  and  Nagamatsu,   and  Arthur,   a  computation  of  the  predicted 

effects  of  theoretical  influence  of  contamination  on  nitrogen  was  made,  using 

values  of  a,  £ ,  and  a/a  of  unity  and  perfect  gas  assumptions  for  nitrogen. 

The  results  are  shown  in  Figure  53,  for  various  size  particles  typical  of 

cluster  sizes  for  CO  ,  indicating  reasonably  good  agreement  with  their  data. 

This  is  taken  to  indicate  some  substantiation  of  the  above  calculations  for 

contamination  in  C0„,  perhaps  within  one  or  two  orders  of  magnitude  of  the 

calculated  contamination  levels  needed  to  influence  homogeneous  condensation. 
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V.    SUMMARY  OF  FINDINGS  AND  CONCLUSIONS 

Ao    Experimental  Findings 
lo  Consistent  data  regarding  the  supersaturation  of  CO  at  onset  of  conden- 
sation plus  the  shape  of  pressure  profiles  has  been  obtained  for  three 
different  nozzles. 

2.  Interferogram  photographs  and  an  approximate  theoretical  estimation 
show  agreement  that  the  boundary  layer  in  the  nozzle  throat  was  of  the 
order  .006"  and  substantiates  the  treatment  of  flow  in  the  mid-stream  as 
1-dimensional  isentropic  flow. 

3.  Departures  from  perfect  gas  behavior  by  CO   in  the  P-T  range  tested  are 
significant,  and  must  be  accounted  for  in  the  calculation  of  local  conditions 
along  an  isentropic  expansion.   Density  and  static  pressure  measurements 
plus  total  pressure  measurements  for  non-condensing  flow  indicate  that  the 
Plank  Equation  of  State  gives  a  good  description  of  low  temperature  C0_ 
behavior. 

U.  Rates  of  CO  expansion,  as  determined  by  nozzle  geometry  have  a  large 
influence  upon  the  degree  of  supersaturation  obtained  prior  to  occurrence 
of  condensation. 

Bo    Application  of  Nucleation  and  Condensation  Theory 

1.  Classical  nucleation  and  drop  growth  theory,  as  applied  herein,  will 
reasonably  predict  the  conditions  at  which  onset  of  CO   condensation  occurs, 
if  flat  film  liquid  surface  energies  are  corrected  by  amounts  within 
suggested  ranges  of  uncertainty  for  small  droplets  or  crystals.   In  these 
calculations  a  droplet  density  equal  to  that  of  CO  liquid  at  the  triple 
point,  1.18  gm/cm  ,  has  been  assumed. 

2.  The  appropriate  correction  to  surface  energy  is  dependent  upon  the 
values  of  mass  and  thermal  accommodation  coefficients,  £  and  a,  and  has  a 
correction  uncertainty  of  about  21$  of  flat  film  values  for  maximum  possible 
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uncertainty  in  £  and  a0   (-  =  1.1T  for  £  =  a  =  0.0:  -  =1.38  for 

a  o 
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E,  =  a   =  1.0.)   Consequently,  the  applicable  value  of  the  parameter 
(o3/pT2)  to  be  applied  to  the  nucleation  rate  equation  corresponds  to 
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values  of  (o 3/p2 ) (pT 2/o 3)  which  lie  between  the  limits  (l.l?)3  and 
(l.38)3where  p   =  1.18  gra./cm.3  and  o   is  the  linear  temperature  depen- 
dent  function  c^  =  5.00  -  ,193T(°C)  dyne/cm. 

3.  Pressure  profiles  obtained  experimentally  indicate  that  considerably 
less  moisture  is  forming  than  that  necessary  to  return  the  mixture  of 
vapor  and  condensate  to  near  equilibrium  values,  or  zero  supersaturation 
downstream  of  the  onset  of  condensation.   Values  of  £  and  a  necessary  to 
produce  agreement  between  the  theory  and  experiment  in  this  respect  require 
that  one  or  both  of  these  parameters  be  considerably  less  than  what  are 
thought  to  be  minimal  values  (about  .01,  for  each.) 

U.  Even  at  their  best  empirical  fit,  classical  nucleation  and  drop  growth 
theory  as  applied  herein  do  not  seem  to  provide  any  means  of  describing 
the  exact  shape  of  the  pressure  profiles  obtained  experimentally  or  the 
extent  of  effects  of  geometry  noted  experimentally,   These  failings,  in 
addition  to  that  of  B3,  above,  are  thought  to  indicate  inadequate  or  incorrect 
assumptions  in  the  treatment  of  the  drop  growth  calculations. 

5.  Application  of  the  so-called  "gassification"  corrections  to  classical 
nucleation  theory  results  in  corrections  in  the  wrong  direction,  unless 
even  larger  corrections  in  the  surface  energy  of  formation  of  CO  ice 
crystals  (than  indicated  above)  can  be  established.   Use  of  solid,  rather 
than  liquid,  droplet  density  in  the  parameter  (o3/p2)  further  increases 
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the  necessary  corrections  to  surface  energies  required  to  provide  agreement 
with  experiment. 

6.  The  influence  of  dirt  particles  or  other  foreign  particles  on  the 
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condensation  of  CO  ,  and  probably  most  other  gases,  may  be  completely 
disregarded  in  supersonic  nozzles,  unless  the  particles  result  from  the 
precondensation  of  another  gas  in  the  stream.   Amounts  of  micron-sized 
particles  necessary  to  influence  CO   condensation  are  about  five  orders 
of  magnitude  in  excess  of  upper  atmospheric  limits  of  concentration,  and 


particles  of  the  order  of  20  a  are  necessary  in  excess  of  about  10" 

3 
particles/cm.   to  influence  homogeneous  nucleation  and  condensation, 

Testing  the  condensation  of  nitrogen  with  pre-condensed  particles  of  CO, 
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in  quantities  mentioned  by  Willmarth  and  Nagamatsu,  and  Arthur,  yields 

theoretical  predictions  which  are  in  reasonable  agreement  with  their 
experimental  observations . 

7.  Additional  knowledge  of  thermal  and  mass  accommodation  coefficients 
are  necessary  to  allow  better  testing  of  homogeneous  nucleation  and 
condensation  theory.   In  addition,  comprehensive  testing  of  the  theory 
will  await  some  better  establishment  of  surface  energy  corrections  for 
small  particles.   Until  this  information  becomes  available  and  an  estab- 
lished treatment  of  nucleation  and  drop  growth  can  be  applied  with  con- 
fidence to  many  vapors,  "recipe"  fitting  may  provide  the  best  means  for 
predicting  nozzle  condensation  behavior. 
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VI .   SUGGESTIONS  FOR  FURTHER  STUDY 

One  of  the  greatest  difficulties  in  application  of  experimental 
results  to  a  meaningful  test  of  theoretical  predictions  results 
from  the  uncertainties  in  value  of  many  of  the  important  variables 
and  parameters  contained  in  nucleation  and  drop  growth  theory.   Con- 
sequently, any  methods  or  devices  of  experimentation  which  would 
decrease  any  of  the  uncertainties,  or  the  effects  of  uncertainty, 
of  such  quantities  as  o,  p  ,  a,  £,  local  temperature  and  the  gas 
dynamics  in  general  would  be  of  great  assistance.   One  means  of 
proceeding  in  this  direction  is  to  contain  the  test  gas  within  large 
quantities  of  another  gas  which  is  non-condensing  in  the  region  of 
interest  of  the  test  gas,  and  whose  gas  properties  are  well-known, 
preferably  being  as  near  those  of  a  perfect  gas  as  possible.   An 
immediate  benefit  of  this  is  to  provide  a  means  of  accurately  determin- 
ing local  properties  at  all  points  of  the  isentrope  with  ease,  since, 
the  temperature  of  the  test  gas  should  be  that  of  the  ideal  "carrier" 
gas,  which  dominates  the  gas  dynamics.   Second,  the  influence  of 
thermal  accommodation  is  greatly  reduced  due  to  the  fact  that  the 
huge  preponderance  of  incident  molecules  upon  the  surface  of  a  formed 
droplet  are  those  of  "carrier"  gas,  and  hence  are  reflected.   The  drop 
is  consequently  kept  at  or  nearly  at  environmental  temperature 
regardless  of  the  value  of  a.   Hence,  there  are  advantages  to  be 
gained  by  testing  gases  for  which  these  uncertainties  exist  (the 
huge  preponderance  of  all  gases)  in  a  carrier  stream  of  such  gases 
as  nitrogen,  helium,  or  dry  air.   Several  complications  arise  from 
this  method,  including  the  need  for  very  precise  measurement  of  mass 
fraction  of  test  gas  in  the  mixture,  and  more  sensitive  instrumentation 
to  detect  the  lessened  effects  of  condensation  upon  the  gas  properties . 
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Highly  steady  flow  systems  would  aid  in  obtaining  the  accuracy 
necessary . 

With  regard  to  some  of  the  uncertainties  pertaining  to  CO  ,  it 
would  be  of  considerable  interest  to  determine  whether  or  not  the 
experimental  behavior  of  condensation  is  influenced  by  establish- 
ing liquid  rather  than  ice  formation  of  condensate  particles.   This 
could  presumably  be  accomplished  by  raising  pressure  levels  to  a  point 
where  condensation  occurs  above  the  ice  point,  75.1  psia.   This 
would  necessitate  stagnation  pressure  in;the  range  of  175  -  250  psia, 
mass  rates,  of  flow  and  a  source  of  00  which  ar^  proportionally 
larger,  heavier  construction  of  some  important  fixtures  in  the  system, 
(such  as  the  stagnation  tank)  and  a  means  of  pressure  profile  procure- 
ment other  than  manometer  board. 

It  is  felt  that  interferometric  means  of  study  offer  enough 
advantages,  particularly  with  more  toxic  or  reactive  gases,  to  warrant 
further  application  to  study  of  condensation  in  nozzles.   The  inability 
to  distinguish  interference  lines  in  the  upstream  region  of  the  nozzle 
encountered  at  high  pressure  levels  in  this  investigation  could 
possibly  be  solved  by  facilitating  higher  camera  shutter  speeds 
(1/150  sec.,  the  maximum  possible,  was  used  in  these  studies)  or  by 
establishing  wall  temperatures  which  are  close  to  stream  stagnation 
temperature.   Two-dimensional  nozzles,  where  interferometric  study 
could  be  successfully  employed,  would  decrease  the  problems  of 
instrumentation  and  would  be  readily  adaptable  to  studies  of  light 
scattering  and  polarization,  yielding  experimental  information 
regarding  droplet  sizes.   Interferometric  photos  involving  conden- 
sation, although  not  applied  in  this  report,  were  obtained,  and  showed 
interference  bands  very  clearly  within  the  region  of  condensation, 
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this  region  being  marked  plainly  by  the  decreased  light  intensity 
resulting  from  the  presence  of  the  ice  cloud. 

Light  scattering  or  other  techniques  which  yield  infor- 
mation regarding  the  average  size  of  the  ice  particles  at  their 
initial  appearance  and  after  subsequent  downstream  growth  would  be 
very  useful.   Since,  as  previously  mentioned,  there  is  some  indication 
that  nueleation  is  the  major  contributor  to  moisture  in  its  early 
stages  of  formation  for  C0o,  this  information  might  be  directly 
related  to  the  critical  cluster  size  of  CO  . 

Additional  investigation  of  crystal  thermal  and  mass  accommoda- 
tion coefficients  would  be  of  help  in  further  establishing  the 
apparent  indication  from  this  investigation  that  either  one  or  both 
have  very  low  values  for  crystalline  CO,,  at  low  temperatures  .   The 
carrier  gas  technique  previously  mentioned  could  isolate  the  influence 
of  a  and  might  facilitate  a  better  means  of  studying  £  and  its  role 
in  condensation. 
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Condensation  and  Experimental  Data,  Nozzle  II  (2.20°) 
Pressure  Tap  Plate  No.  1.  a  =  1.00,  o/o^   =  1.18, 
£  =  .00008 

Figure  U8   Comparison  of  Empirically  Fitted  Theoretical  A66 

Condensation  and  Experimental  Data  -  Nozzle  II  (2.20°) 
a  =  1.0,  C  =  .00008,  0/0     =  1.18 

OO 

Figure  k9     Comparison  of  Empirically  Fitted  Theoretical  A67 

Condensation  and  Experimental  Data  -  Nozzle  II  (0.80°) 
a  =  1.0  £  =  .00008,  0/0     -   1.18 

OO 

Figure  50  Theoretical  Influence  of  Contamination  on  Condensation    A68 
of  CO 

Figure  51  Theoretical  Influence  of  Macroscopie  Contamination  on     A69 
Condensation  of  CO  ,  Using  Theory  Empirically  Fitted  to 
Homogeneous  Condensation  Data 

Figure  52  Theoretical  Influence  of  Microscopic  Contamination  on     A70 
Condensation  of  CO2,  Using  Theory  Empirically  Fitted  to 
Homogeneous  Condensation  Data.  £  =  .00008,  0/0     =1.18, 
pT  =  1.18  GM./CM 

Li 


Aiv 


Figure  53  Theoretical  Effect  of  CO  Contamination  on  Nitrogen       ATI 
Condensation.   Supersaturation  at  .1%   Moisture 
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FIGURE    7,    ASSEMBLY   VIEW  AND  GEOMETRY  OF   NOZZLE    I, 
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FIGURE  9,  NOZZLE  I;  DISASSEMBLED,  SHOWING  OPTICAL  GLASS  AND  METAL 
PRESSURE  TAP  PLATE  SIDE  WALLS,  AND  METAL  CLAMP  BARS. 
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FIGURE  10,  NOZZLE  II;  DISASSEMBLED,  SHOWING  HINGED  UPPER  AND  LOWER 
WALLS,  AND  PRESSURE  TAP  PLATES  1  AND  2 
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(A)  NO  FLOW 


(B)  SCALING  PLATE  (1.00") 


(C)  FLOW,  SHOWING  SHIFT  IN  INTERFERENCE  BANDS 
FIGURE  II,  INTERFEROMETRIC  PHOTOGRAPHS  FOR  TEST  1  OF  DENSITY  PROFILES 
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(A)  NO  FLOW 


(B)  SCALING  PLATE  (1.00") 


(C)  FLOW,  SHOWING  SHIFT  IN  INTERFERENCE  BANDS 
FIGURE  12,  INTERFEROMETRIC  PHOTOGRAPHS  FOR  TEST  2  OF  DENSITY  PROFILES 
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FIGURE  17,  EXPERIMENTAL  NORMAL  SHOCK  CHARACTERISTICS  OF  C02 
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FIGURE  18  ,  PREDICTED  NORMAL  SHOCK  CHARACTERISTICS  OF  C09  -  PLANK  EQUATION 
OF  STATE. 
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FIGURE  19,  COMPARISON  OF  PREDICTED  C02    ISENTROPES 
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FIGURE  20  ,  COMPARISON  OF  PREDICTED  C02  ISENTROPES 
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FIGURE  21  ,  CONSTANT  NUCLEATION  RATE  PROFILES  OF  CLASSICAL 
NUCLEATION  THEORY  FOR  CO„ .  %      =1.00. 
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FIGURE  22,    EXPERIMENTAL  PRESSURE   RATIO   PROFILES.    NOZZLE   I, 
PRESSURE   PLATES    1   &  2 
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FIGURE  31,  CONDENSATION  DATA. 

LOCAL  CONDITIONS  AT  ONSET 
OF  CONDENSATION. 
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FIGURE  36,  ONSET  OF  CONDENSATION,  NOZZLE  II, 

COMPARISON  OF  THEORY  AND  EXPERIMENT 
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FIGURE  37,  INFLUENCE  UPON  THEORY  OF  USE  OF  LIQUID 

VERSUS  SOLID  DROPLET  DENSITY.  NOZZLE  11(2.20°) 
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FIGURE  42,  THEORETICAL  OCCURRENCE  OF  .  170  MOISTURE  FOR  a  =  £  =  1.00, 
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FIGURE  43,  THEORETICAL  OCCURRENCE  OF  .1%  MOISTURE  FOR  a=  £  =  .01, 
CORRECTED  SURFACE  TENSION  vl 'a_  =  1.20 
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APPENDIX  B 

REVIEW  OF  NUCLEATION  RATE  EQUATION  DEVELOPMENT 

The  following  is  a  review  of  the  classical  nucleation  rate  equation  as 

(19) 
developed  by  J.  Frenkel    .  Although  Frenkel  was  not  the  first  to  investi- 
gate nucleation  theory,  his  result  is  not  greatly  different  from  earlier 
authors  mentioned  in  the  text,  and  his  method  provides  relatively  straight- 
forward consideration  of  the  important  aspects  and  assumptions  in  classical 
nucleation  theory 

Frenkel *s  droplet  model  was  based  on  Volmer's    earlier  model  of  drop- 
lets as  generalized  molecules,  neglecting  intermolecular  forces,  and  his 
method  of  development  was  based  on  that  of  Zeldovich    .  He  showed  that 
the  distribution  of  clusters  of  vapor  molecules,  resulting  from  the  kinetic 
interaction  of  the  molecules  could  be  described  by  a  Boltzmann  distribution. 
That  is,  that  (in  equilibrium) 

N(g)  =  No  exp  (-M/KT)  (B-l) 

where     N  =  total  number  of  molecules 
o 

$   =  H  -  TS 

A<J>  =  $T  -  $  =  change  in  free  energy  of  a  group  of  (g)  free  mole- 
cules  when  forming  a  liquid  embryo. 

Several  simplifying  assumptions  of  debatable  validity  are  made  at  this 
point.   First,  the  droplet  model  is 'taken  to  be  a  liquid  sphere,  having  sur- 
face energy  described  by  the  flat-film  liquid  surface  tension  (or  some  cor- 
rected function  thereof) .   Second,  the  total  change  in  free  energy  is  assumed 
to  be  the  work  of  formation  of  the  surface  of  the  droplet,  plus  the  change 
in  free  energy  resulting  from  the  change  in  state  of  g  molecules  condensing 
from  vapor  to  liquid,  with  the  spherical  droplet  assumed  to  be  stationary 
in  the  vapor  system,  and  having  no  rotational  or  vibrational  energy.  This 
is  the  important  distinction  between  classical  and  more  recent  "gassified" 
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droplet  models  of  nucleation  theory.   Third,  the  droplet  is  assumed  to  form 
at  environmental  temperature,  an  assumption  which  may  also  be  invalid    . 

Using  these  assumptions,  and  defining  4>  as  the  free  energy  per  molecule, 
A$  can  then  be  writeen: 

A*   =   g(4>L   -   <J>V)    +   W2o  (B-2) 

=   -  T(AS)    +   W2o 

=  -  gKT  £n  ^~  +  h-nr2o  (B-3) 

where  p^  equals  the  flat  surface  equilibrium  pressure  of  a  formed  embryo  of 
(g)  molecules  and  p  is  the  existing  local  pressure. 
Hence, 

N(g)  =  NQ  exp(g  £n  £  -  ^|~)  (B-U) 

Saying  that: 

gm  =  ^^  PL  (B-5) 

N(g)  -  NQ  exp  {  g  to  J  -  j£2.  (^)2/3}  (B-6) 

00 

=  NQ  exp  {  g  *n  E-  -  (j^)  g2/3} 

00 

u  i,   /3m  N  2  /  3 

where     e  =  4tto(i )   ° 

imp 

Inspection  of  Eq.  (B-6)  shows  that  a  reasonable  equilibrium  distrib- 
ution results  when  p  is  less  than  p  ,  that  is,  in  a  superheated  vapor.   How- 
ever, in  a  supersaturated  vapor,  an  unreasonable  and,  in  fact,  impossible 
equilibrium  distribution  is  predicted  as  seen  in  Fig.  (B-l),  since  it  pre- 
dicts N(g)  ->  oo  as  g  ->  °°  .   Acknowledging  that  an  equilibrium  distribution 
of  drop  sizes  is  unreasonable  under  supersaturated  conditions,  early  authors 
treated  the  system  as  one  in  a  state  of  quasi -equilibrium,  where  droplets 
are  imagined  to  be  removed  after  reaching  some  size  larger  than  the  critical 
radius,  r*,  which  occurs  at  the  maximum  of  the  supersaturated  free  energy 
change.   This  maximum  in  A*  (minimum  in  N(g))  occurs  when; 


yields 
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BU. 


1/3 


-r, 


(A*)   =  0  =  -  KT  In  p~    +  f-  eg     ^ 


D  ?  -1/3 


linen 


Hence 


max        3  3 


A$  =  _  gKT   in  £     +   Uiir2cr  =  -  Trr2a  (B-7) 

max  p  3 


KT(|  *r3  -L)  In  *  -  f  *r2c 
3     m      Poo   3 


in^  = 


2mo 


Poo   PLKTr*  (B-8) 

Equation  (B-8)  is  the  Kelvin-Helmholtz  expression  describing  the  conditions 
under  which  a  liquid  drop  may  remain  in  equilibrium  with  its  vapor  phase. 
The  radius  (r*)  and  the  number  (g*)  of  atoms  of  a  critical  embryo  are  thus 
determined.  A  variation  of  r  (or  g)  in  either  direction  results  in  a  de- 
crease of  free  energy  from  the  maximum,  and  consequently  the  drop  is  in 
unstable  equilibrium  at  this  point . 

In  considering  the  quasi-equilibrium  state,  and  ignoring  interaction 
between  nuclei,  it  is  said  that  nucleation  rate  is  the  net  rate  at  which 
particles  of  a  given  size  (g)  are  formed.  To  maintain  quasi-equilibrium, 
it  is  further  assumed  the  net  rate  at  which  (g+l)- sized  droplets  are  formed 
from  (g)-sized  droplets  is  equal  to  the  rate  at  which (g)- sized  droplets  are 
formed  from  (g-l)-sized  droplets,  and  so  on.   That  is,  nucleation  rate  is 
independent  of  g,  in  order  that  the  number  of  droplets  of  any  size  remains 
constant  and  the  rate  at  which  droplets  of  larger  size  c  be 

removed  is  also  equal  to  the  nucleation  rate.   Correspondingly,  it  is  im- 
agined that  single  molecules  are  added  to  the  system  at  a  rate  equal  to  the 
depletion  due  to  nucleation. 

Using  this  model  then,  and  assuming  that  growth  of  nuclei  is  due  only 
to  interaction  with  free  molecules,  the  nucleation  rate  is 
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J(g  +  1,  t)  =  S  s(g)x  n(g,  t)  -  Y(g  +  1)  s(g  +  1)  n(g  +  1,  t) 

(B-9) 

where  6  =  rate  at  which  molecules  strike  unit  surface  area. 

y(g)   =  rate  at  which  molecules  leave  the  surface  per  unit  area 

n(g»t)  =  distribution  of  non-equilibrium  drops 

s(g)   =  surface  area  of  a  g-sized  drop. 

Assuming  that  no  effects  of  droplet  curvature  exist,  and  neglecting  the 
cluster  motion, 

6  =  (2TrmKT)i/^ 
Also,  it  is  said  that  if  nucleation  rate  were  zero,  then  the  non- 
equilibrium  (or  quasi-equilibrium)  distribution  would  in  that  case  have  to 
be  equal  to  the  equilibrium  distribution  N(g).  This  is  easy  enough  to  see 
in  the  superheated  case,  but  perhaps  a  bit  difficult  to  accept  when  N(g)  is 
the  fictitious  mathematical  result  corresponding  to  supers at urat ion. 
Applying  this,  when  J(g,  t)  =  0 

y(g  +  1)  s(g  +  1)  n(g  +  1)  =  6  s(g)  n(g) 
when      n(g)  =  N(g) 
Hence, 

y(g  +   1)    s(g  +  1)    N(g  +  1)    =   6   s(g)   N(g) 
and  Equation   (B-9)   becomes 

J(«.   t)   =   B  .(g)   N(g)    (  Sifoii  -  Si^-f  )  (B-10) 

Here   some  rather  notable   simplifications   and  manipulations   facilitate 
a  solution  to   (B-10) .      For  large  values  of  (g)   where  these   functions  may  be 
treated   as   continuous,   this   is  approximately: 

J(g,  t)  =  $  s(g)  N(g)  ^  (£) 
Hence,  since  it  is  assumed  that  J(g,  t)  =  J(t) 

/^■-f/^ite  (B-11) 

Using  the  boundary  conditions 
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(a)     Jo  i.o  at  g  =  1  (b)      Aim     J  =   0 

g      -*      00 


0  ne   obtains 

j  ■  far  <B-i2> 

J  !    sN 

ll  r  a<j>  I 

Since  —  =  —  exp    ^  —  J    has   a  sharp  maximum  in  the  vicinity  of  g  =  g* ,    it 

o 

is  possible  to  approximate  (B-12)  as: 

6  s(g*)  N 

J  =  — 2 

/Texp  (A*/KT)  dg  (B-13) 

Using  Equations  (B-5),  (B-6)  and  (B-8),  the  exponent  may  be  written: 

M  =  _  g  £n  £_  +  e   2/3 
KT     S  *p    KT  g 

=  _  -    r       2mo     I   +  itH2.  r  3m     -[2/3     2/3 

6    L  pTKTr*  J        KT      Uttpt   j  e 

Li  L 

=  A    [  3(g/g*)2/3   -  2(g/g*)'] 

where 


WpL\l/3     *2/3    Wr*2 
p  KT    \     3m    /         g  3KT 


am      fZlLV 
™   \     3m    / 


ft   s(g*Q    No 


J  =   /"  exp    (  A    [  3(g/g*)2/3  -   2(g/g*)]   }      dg~~  (B-lU) 

This    integral  occurs   in  the  Becker  and  Doring,    and  Frenkel   derivations, 
and  can  be  solved  approximately  with  the   introduction  of 

0    =   (g/g*)1/3-  1  (B-15) 

A 
Since  the  integrand  of   (B-lM    has   a  very  sharp  maximum  e     at   g  =  g* ,   the 

important  values  of  ft  are  small,    and  the   limits   of  the   integral  may  be  taken 

at  +°°  and  -°°. 

Also,    since  (g/g*)2/3  =  ft2  +   2ft  +   1 

dg  =   3(1  +  ft)2  g*dft 

And  3(g/g*)2/3   _   2(g/g*)    =   -   3ft2   -   2ft  +   1 

The   integral  may  be  written  approximately 

3g*eA   j7ro   (1  +   ft)2  exp   (-   3Aft2   -  2Aft3)    dft 

-   3g*eA    r     exc   (-   3Aft2)    dft 

'  —  oo 


BT. 


3g*e  (tt/3A)  (B-16) 


Hence,  the  nucleation  rate  finally  becomes 
* 


— T*~°L  "57  J      exp(-A)  (B-17) 


g        L   3tt 


with 


N     =  p/KT 


6      =   p/(2TrmKT)l/2 
s(g   )   =     -Q- 


and.  with  the  definitions   of  e   and  A,    (BIT)    can  be   shown  equivalent   to 

j  =   (E_f  5L  (gg.)l/2exp(-U™r*2) 

4t;    p    W        exp*      3KT      ;  (B_lQ) 

This  is  an  expression  of  the  number  of  nuclei  formed  per  unit  volume 

per  unit  time,  due  to  the  existence  of  a  non-equilibrium  distribution  n(g). 

By  assuming  (j)  was  independent  of  (g)  as  Becker  and  Doring  did,  it  was 

evaluated  using  only  the  boundary  conditions  of  n(g)  and  the  properties  of 

A$(g). 

Discussion  and  criticism  of  Eq.  (B-17),  and  of  the  importance  of  its 
various  constituent  values  is  contained  in  the  text. 
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APPENDIX  C 

REVIEW  OF  DEVELOPMENT  OF  DROP  GROWTH  EQUATIONS 

The  following  development  follows  the  method  of  Oswatitch    in  consid- 
ering the  necessary  conservation  of  mass  and  energy  of  growing  droplets. 
Included  are  a  correction  of  average  incident  energy  to  be  2RT  rather  than 
3/2RT,  and  with  generalized  non-unity  values  of  both  mass  and  thermal  accom- 
modation. An  important  assumption  is  that  molecular  velocity  distribution 
are  Maxwellian  before  and   after  collision  and  accommodation  with  the  drop- 
let surface. 

Consider  a  drop  of  temperature  T   in  an  atmosphere  p  and  T.   It  has  an 
incident  mass  flux  equal  to 

P  ,  .. 

(2ttRT)1/z 

with  a  proportion  £  (condensation  coefficient)  condensing  and  the  rest  being 

reflected.  Since  the  drop  is  very  small,  its  temperature  is  assumed  to  be 

constant  throughout.   The  mass  transfer  from  the  drop  is  equal  to  that  which 

would  be  incident  in  an  environment  T  and  p  ,  the  equilibrium  values  of 


temperature  and  pressure,  where 
Jin 


PD  2ma 


Hence,   the  mass   flux  per  unit  area  passing  out  through  the  surface   is 

(2irRT    W2 
Thus,   the   conservation  of  mass  may  be  written 

*Lkvr2    ( dT}    =   U7ir2  [  (2  J?)  1/*    "    (2irRT°)l/*] 

or  dr  -      €  (  P.  ...  PD      \  /  n   ,  \ 

dt   -(2irR)l/2    I    tT7T-^I72J  ^'^ 

assuming  that   C    is   constant,    or  at  most,   a  function  of  the  drop  temperature 
only. 

Considerations   of  energy  conservation   are  made  after  first   defining  the 
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the  thermal  accommodation  coefficient,  involving  the  change  in  temperature 
of  reflected  incident  molecules. 
T   ...   -  T 


a 


"refl.  (Thermal  Accommodation 


T  -  T  Coefficient) 


(C-2) 


The  energy  flux  to  the  drop  is,  per  unit  area 

(a,m)Vii  [  2RT  ] 

and  the  energy  leaving,  upon  reflection 

7^%  {*  [  T  ♦  a(TD  -  T)]J 
In  adaiiion,  the  evaporating  mass  flux  carries  energy 

(2ttRTd)V^  f  2RTD  ] 

assuming  the  evaporating  mass  flux  also  to  be  Maxwellian. 
The  rate  of  change  of  droplet   internal  energy  is 
h        r>  D    |   o  dr 

3  *r  V  3T  +  k"    at  pl  ufD 

UfD  '  hfD  "  W  RT  "  "fg 
Hence,  conservation  of  energy  requires, 
)        dTn        a 

1  "  V  3T  +  h"2  f  pL(7^T  RT  -  V  ■ 

B(2RT)    -   C3D(2RTD)    -   {l  -   O    32R    [  T  +   a(TD   -  T)] 
or  rp    c      dT  6n  Tn  T 

-5?&  [7^TRT-hfJ 

"here  8  =   (2^l/2  BD  "  ^uKr")'/* 

An  order  of  magnitude  analysis   shows  that  the  left  hand  term  is   small 

compared  to  the   last   term  on  the   right   hand  side.      Taking  typical  values; 

r  ~  10~7  cm  h  JC  T  ~  3 

fg     P 

c  ~  .3  cal/gm  °K         R  ~.0U5  cal/gm  °K       T  ~  250°K 
Y/(y  -  1)~);  dT  /dr  —  7  x  107  °K/cm 


C3 


dT 
re   D 


dr 

=  .008  ~10 


hf^ 
Y  -  1   C 


This  signifies  that  the  energy  required  to  change  the  droplet  temper- 
ature during  the  expansion  is  much  less  than  the  latent  heat  imparted  from 
condensing  incident  molecules.   Neglecting  this  term,  the  energy  equation 
becomes  ? 

o  =  2€(i .  JaS>> .  2d  -  c)  .(Ji  - 1)  ♦  Stte  (_jl_,(A  .  d 


where 

C  T 
P 

Since 

dr        E      /„ 

dt       pL 

•  V 

(c-i: 


The  energy  equation  can  he  further  reduced  to 

o  =  (i  -  |2  lD-  )  -  „(i  -  i){^  -  D  +  |(i  -  ^H^fyH*  -  i) 

(0-3) 

Solution  to  the  drop  growth  equations  involves  an  iterative  determin- 
ation of  drop  temperature,  T  ,  from  Eq.  (C-3)  and  the  resulting  value  of 
growth  rate  from  Eq.  (C-l). 

An  interesting  indication  of  the  importance  of  considering  both  mass 
and  energy  balance  may  be  seen  by  closer  inspection  of  Eqs .  (C-l)  and  (C-3) 
It  is  seen  that  the  requirement  for  a  drop  to  grow  is  that 

££ .  M  (i  -  %  >  i        ...  £  «  1  (C-U) 

dt   p       6  (3 

Li 

Re-arranging  Eq.  (C-3)  into  the  form: 


3D 


[^r  +  l(7^T)(x  -1'  ]  =  i  -a(i-D(/-i)  *|(7fi)(»-i) 


we  obt  ai  n 


cu. 

It  is  clear  that  whenever  T   >  T  the  denominator  of  Eq.  (C-5)  is  always 
greater  than  the  numerator .   Hence  from  Eq.  (C-H)  the  droplet  will  grow. 
Even  in  the  most  unfavorable  case,  that  of  a  =  0,  where  the  drop  temperature 
becomes  much  hotter  than  the  surrounding  environment,  the  value  of  8^/R  will 
remain  less  than  unity  and  the  drop  will  grow.   This  is  really  a  nroof  of 
the  intuitive  fact  that  a  droplet  of  size  greater  than  critical  radius  must 
continue  to  grow  in  a  supersaturated  environment.   Conversely  when  T   <  T, 
3/3  is  greater  than  unity  and  the  drop  evaporates,  supnortinp;  our  knowledge 
that  evaporating  liquids  are  cooler  than  their  superheated  environments. 
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APPENDIX  D. 

MODIFICATION  OF  FOREGOING  GAS  DYNAMICS,  NUCLEATION  AND  DROP  GROWTH 
EQUATIONS  FOR  DEPARTURE  FROM  PERFECT  GAS  DURING  AN  ISENTROPIC 
EXPANSION. 

The  following  corrections  introduced  into  the  perfect  gas  relations  to 
account  for  departure  from  perfect  gas  behavior  are  based  on  the  empirical 
observation  that  compressibility  factor  z  =  Pv/RT  is  almost  constant  during 
an  isentropic  expansion.   This  is  illustrated  by  Table  15  which  lists  the 
values  of  P,  T  and  z  for  all  of  the  condensation  tests  of  this  investigation, 
Changes  in  z  generally  encountered  were  less  than  .5%   from  stagnation  press- 
ure to  maximum  supersaturation  prior  to  condensation.  This  in  effect  allows 
introduction  of  "effective"  perfect  gas  relations, 

Pv  =  R'  T  where  R'  =  zR  (D-l 

and  the  gas  dynamics  are  affected  accordingly: 

Tds  =   0  =  C  dT  -  vdp 
P 


n     pdv   vdp     , 

_  _L_  [  dP    dpi   Y-l  dp 
Y-l  L  P     p  J    Y   p 


Hence     dp_     dp_ 
P      p 


l-z(^)  .   *>(!,)  (D_2 


p    y 

Hence,  the  isentropic  relations  of  a  perfect  gas  may  be  modified  by  replacing 
Y  with  Y*  as  defined  in  equation  (D-2). 

Similarly,  the  free  energy  of  formation  of  a  nucleus  is  modified  as 
follows.   The  free  energy  change  from  vapor  to  liquid  of  (g)  molecules  (added 
to  the  free  energy  of  surface  formation  in  the  classical  droplet  model)  is 
d<t>  =   dH  -  d(TS) 
If  it  is  assumed  that  the  nucleus  forrs  at  environmental  temperature. 
dH  =  SdT  =   0 


d4>   =   -TdS   =   vdp  =   zRT 

P 


dp_ 


D2. 

4>2   _   <|,  x   =   zgkT   *n  —  =  -   zgkT   An  ^ 
^1 


P,  Poo 


or,  the  total  change  in  free  energy  of  formation  of  a  classical  model  droplet 
becomes 

$2   -  *i   =  -zgkT  ^n  £  +  Uirr2  o 

Poo 

By  proceeding  as  in  Appendix  B,  it  can  be  further  shown  that 

A<I>     ■   -J  iir2  0   at   r  =  r* 
max      3 


and 


An  £  = 


2m° 


Poo    PL  zKTr* 
These,  plus  the  corrections  to  the  drop  growth  equations,  may  be 

introduced  to  the  perfect  gas  equations  by  the  simple  substitutions 

i   m 
m  "  z  (D-3) 

y1   "  1  +  (Y  I  i  )(i  .  z)  (D-U) 

In  a  computation  involving  condensation  of  a  "real"  gas,  it  is  necessary 

to  first  determine  the  real  isentrope,  determined  through  use  of  the  equation 

of  state  as  mentioned  in  the  text. 

From  this,  by  use  of  the  equation  of  state  and  Equation  D2,  respectively, 

z  and  Y'  are  established  at  all  points  of  the  expansion  and  it  is  then  a 

simple  matter  to  determine  the  value  of  y   to  be  entered  into  the  computation. 

The  value  of  z  is  taken  as  the  mean  of  stagnation  and  condensation  values 

(indicated  in  Table  lU)  and  a  linear,  temperature  dependent  approximation  to 

Y  between  the  same  two  points  is  made, 

Y  =  Y   -  Y.  T(°C) 
a    b 

where  values  of  Y   and  Y,  are  given  in  Table  15. 
a      b 

Consequently,  the  computation  scheme  treats  a  value  of  y  which  will  yield 
the  real  isentropic  exponent  Y'  when  treated  with  the  corrections  (D-3)  and 
(D-U).   This  matches  the  isentrope  prior  to  condensation  with  the  best 
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estimate  of  local  conditions  and  calculates  corrected  nucleation  and  drop 
growth  rates  within  the  theoretical  framework  of  perfect  gas. 

Perfect  gas  calculations  may  be  made  by  simply  setting  z  =  1.00,  and 
entering  the  desired  value  of  y   =  y    ,  y     =  0. 

cL     0 

It  might  be  pointed  out  that  the  compressibility  factor  begins  to  vary 
more  rapidly  after  condensation  commences,  and  that  the  computation  scheme 
ignores  this0   Some  error  in  profile  shapes  may  result,  but  probably  are  no 
more  in  doubt  than  those  resulting  from  perfect  gas  treatment.   The  value 
of  this  treatment  lies  in  its  attempt  at  establishment  of  the  best  estimate 
of  conditions  and  behavior  at  the  point  of  onset   of  condensation. 

Another  possibility  for  the  consideration  of  real  gas  behavior  in 
condensation  would  be  to  completely  rewrite  the  entire  processes  of  nucleation, 
drop  growth  and  gas  dynamics  in  terms  of  an  equation  of  state,  as  was  done 
for  the  non-condensing  isentropes  and  total  pressure  computations.   This 
could  only  be  done  by  greatly  increasing  the  complexity  of  an  already  complex 
computation  and  at  the  cost  of  excluding  the  many  perfect  gas  relationships 
which  so  greatly  simplify  the  mathematical  treatment  of  condensation  in  gas 
flow  processes. 
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APPENDIX  E,  A  CLOSED  FORM  SOLUTION  TO  THE  MOMENTUM  INTEGRAL  BOUNDARY  LAYER 
EQUATION  INCLUDING  APPROXIMATE  EFFECTS  OF  STREAMWISE  DENSITY 
VARIATION 

In  most  cases  involving  the  study  of  boundary  layers,  effects  of  changing 
density  are  negligible „   However,  in  such  cases  as  nozzle  flow,  where  sonic 
velocities  are  approached  or  exceeded,  acceleration  forces  due  to  the  presence 
of  the  density  gradient  become  appreciable,  and  near  the  throat  of  a  supersonic 
nozzle  are  of  the  same  order  as  those  due  to  area  change  in  their  influence 
on  main  stream  velocity  and  hence  boundary  layer  growth <>   A  study  of  the 
momentum  integral  equation  has  revealed  that  a  closed  form  solution  is 
possible  for  both  laminar  and  turbulent  flow.   In  the  laminar  case,  an 
extremely  simple  form  of  solution  is  possible,  enabling  one  to  determine 
estimates  of  boundary  layer  growth  and  depression  in  supersonic  nozzles, 
using  only  slide-rule  calculations. 

If  it  is  assumed  that  density  and  pressure  vary  in  the  direction  of  flow, 
but  not  across  the  boundary  layer,  the  momentum  integral  equation  may  be 


written  as  follows : 
6 


dU 


—  p/   (U-u)udy  +  f>  —  /   (U  -  u)dy  =  T 


(E-l) 


">| 

A(x) 

-> 

U(x) 

6 

W 

p(x) 
\ 

-  =  f(£) 

u     V 


*  y 


For  the  general  case,  it  is  assumed  that 
u/U  =   f(y/6) 

Vpu2  ■  a,  iy Vn 


(E-2a) 
(E-2b) 


Then,  Equation  (l)  may  be  reduced  to: 


E2. 


b  §-(p6u2)  +  b  P6U^  =  a,  (^  PU2 
1  dx  2     dx     1   U<5y 


(E-3) 


where  bj,  b2  and  a,  are  dependent  upon  the  assumed  profile. 
Equation  (E-3)  may  be  manipulated  to  appear  in  the  form: 


+   6' 

?   dx 


2b  +  b  JTT   b   jn 
1    2  dU    i  dP 

U    dx    P  dx 


=  a 


1  MT*; 


or 


d(62) 
dx 


62 


£i  dU   2  dp 
U  dx  "  P  dx 


■   C 


2  V 


7n 


(E-U) 


where 


2(2b  +  b  ) 
C  =  1 2_ 


2a 
2     Di 


1        bi 

Equation  (E-U)  may  be  put  into  integrable  form  by  using  the  integration 
factor: 


/Pdx  du     dD 

m  =   C     ;   Pdx  =   C  f^-   +   2^- 
'  1  U      p 


Hence,  m  takes  the  convenient  form: 
C 
m  =  U  1  p2 

C 

ru  ^  l     p 

or„  more  conveniently    m  =  {—   J    r— — "j 


(E-5) 


and  the  new  form  of  the  differential  equation  is 


^  (m62)   =  mQ 


where   Q=  C2  (^-)/n  6 


(E-6) 


Hence 


d7  m   )   "  C2  m(U}  '  6 

n-1 


n-1 
2n '      Q'"   (m«2)  2n 


/ ■,    n-lv     i  / 
,(1  -  -^r)   rv  /n 


,n-l, 

and:   (m*2)"  2n   f-(m62)   =  c.rfh    •  (^)/n 

dx  z        U 


E3« 


Integrating,  one  finally  obtains 


(m62) 


n+1  i 
2v  2n 


C  (n+1) 
2n     J 


il 


x  m  2n   (-)'    dx 
o 


(E-7) 


Equation  (E-7)  is  the  most  general  form  of  solution  to  the  momentum  integral 
equation,  approximately  valid  for  a  compressible  gas  in  either  a  laminar  or 
turbulent  boundary  layer  <, 

It  should  be  noted  that  the  integrating  factor  m  can  often  be  more 
conveniently  expressed  as  follows: 


pdx     = 


c:<-r-  r'  +2 


s 


m       =      e'Pdx     m     A-Cl   p2-Cl 


or,  more  conveniently 
m    = 


(a_)"Ci  (£_)2"Cl 

o 


(E-8) 


Laminar  Boundary  Layer 

Considering  Equation  (E-7)  for  the  case  of  a  laminar  boundary  layer, 


where  n   ~  ^«  we  immediately  have: 


_x 


m< 


=   C2/x  m($  dx 
o         o 


(E-7a) 


Assuming  that  viscosity  U  is  nearly  constant  and  substituting  for  m  with 

Equation  (8)  yields: 

C  V        x  m 

62  =  -L_  /   m(I_)  dx  +  _2  62 


m    x    PU' 
o 


m   o 


m 
_c 

m 


=   C  U  (-  )°1  (£_)Cl"2   /X   (A-)'0^^-)  "  (i-)  dx  +  -2  62 
U2M  VA*;    VP  '       Jx      KA*J  Po     PU 

o        o 

Further,  for  a  nozzle:   mass  flow  rate  =  constant  =  PUA  =  P*U*A* 
Hence: 


S2   = 


Cou      Ci     Ci  -2  x   A   1-Ci  n      2-Ci  )     m  6  2 
P»U*  V    VD  '       x   A        o       !     m 


Equation  (E-9)  is  the  specialized  form  of  Equation  (E-7)  for  a  laminar 
boundary  layer.   Reviewing  Equations  (E-3)  and  (E-1+),  it  is  seen  that  the 
constants  0j  and  C2  vill  depend  on  the  polynomial  used  to  approximate 
the  shape  of  the  boundary  layer  velocity  profile. 

For  flow  in  a  nozzle,  (A/A*)  is  known  and  (p/p  )  can  be  easily 
determined.   Hence,  the  function; 

1-Ci     2-Ci 

F  =  m(I*)  =  (I*}      (|H  (E"10) 

o 

may  be  also  easily  determined  and  plotted.   The  integral  may  be  determined 

by  the  area  under  the  curve  F,  and  at  any  point,  x: 

x         m 

S2  =  K,  m"1  /   F  dx  +  -2.  62 
1      x         m   0 
o 

For  <$o  =  o  at  x  =  0°. 

-1   x 
<$2   =  K  m   /   F  dx  (E-9a) 

IX 

o 

where  Ki  =  Pu* 

Equation  (E-6)  was  applied  to  a  nozzle  of  known  geometry,  using  the 
assumption; 

U      6  V 

for  which  Cj  =  9   C2  =  30 

At  first  glance,  the  assumption  of  a  laminar  boundary  layer  may  seem 

(Rl) 
unjustified,,   However,  Launder     has  shown  that  turbulent  boundary  layers 

under  sufficiently  negative  pressure  gradients  will  revert  to  a  definite 

laminar  forrru   Also,  for  the  conditions  of  this  calculation  the  observed 

boundary  layer  at  the  nozzle  throat  was  such  that  the  boundary  layer  Reynolds 

Number  was  about  2000  at  the  throat.   This  is  less  than  the  transition  value 

and  would  indicate  that  a  reversion  to  a  laminar  form  could  have  taken 

(Rl) 
place    o   Hence,  the  laminar  assumption  is  not  unreasonable. 


E5. 

Calculations  were  made  according  to  Table  I  for  CO  at  stagnation  cond- 
itions of  Qh   psia  and  -10  F,  according  to  the  corresponding  values. 

F    =    (A/A*)-8  (P/P  rT  n      oL   . 

o  p   =  oU  psia 

ro        ^ 

30P  o 

K,    =    T£T7*  for  COo  at  T      -10  F 

1        p*U*  z      o 

Results  of  calculations  are  shown  in  Figure  El. 

At  the  conditions  noted  above  and  for  the  nozzle  in  question,  limited 
interferometer  data  provided  visual  means  of  measuring  the  boundary  layer 
thickness,  and  four  values  are  shown  in  Figure  El.   These  observations  are 
estimated  to  be  accurate  ± 15%. 

Also  shown  in  Figure  1,  to  assist  in  estimating  the  relative  size  of 
the  boundary  layer,  are  those  values  of  uniform  flow  flat  plate  laminar 
boundary  layers  calculated  for  the  densities  and  velocities  at  the  entrance 
and  throat  of  the  nozzle,  respectively. 

It  is  seen  that  the  limited  experimental  data  seems  to  agree  well  with 
the  approximate  calculation  of  boundary  layer  thickness.  More  experimental 
data  is  necessary  for  a  more  thorough  evaluation. 

The  extremely  small  boundary  layers  resulting  in  the  pressure  gradients 
which  occur  when  compressibility  becomes  important  make  it  appear  that 
further  interferometer  observations  would  be  an  advisable  method  of  gaining 
such  data  if  it  is  not  readily  available  in  the  literature. 

The  fact  that  the  past  "history"  of  the  boundary  layer  has  a  greatly 

decreased  effect  on  the  value  of  boundary  layer  thickness  in  the  region  of 

rapi^  acceleration  is  quite  apparent  both  from  the  curve  in  Figure  1  and 

from  inspection  of  the  equations.   Equation  (E-9)  shows  that  the  value  of 

62,  which  exists  at  the  beginning  of  integration  must  be  multiplied  by  the 
o 

ratio  m   /m  .   Hence,  if  the  value  of  m  increases  appreciably,  as  it  does 
xo  x 

when  flow  area  A  and  density  p  decrease,  then  the  term  (m  /m  )62  contributes 

xo  x  o 


E7, 

little  when  added  to  the  integral  term. 

The  fact  that  one  is  limited  to  an  assumed  boundary  layer  velocity 
profile  during  an  integration  does  not  restrict  the  possibility  of  breaking 
the  flow  into  regions,  and  applying  Equation  (E-9)  to  different  velocity 
profiles  (i.e.  different  constants  C^  and  C2)  over  the  various  intervals  of  x, 
with  relative  ease. 

Turbulent  Boundary  Layer 

It  is  apparent  that  Equation  (E-7)  is  readily  applicable  to  turbulent 
boundary  layers  in  its  present  form,  merely  by  selecting  the  expressions 
desired  for  Equations  (E-2a)  and  (E-2b),  that  is,  the  velocity  profile  and 
wall  shear  dependence „ 

However,  since  most  persons  are  accustomed  to  referring  to  turbulent 
boundary  layer  behavior  in  terms  of  momentum  thickness  9,  and  shape  factor 
H  =  6*/9,  it  is  worthwhile  to  re-derive  Equation  (E-7)  in  terms  of  6  and  H, 

6  6 

^-  p  /   (U  -  u)udy  +  P  ~     j    (U  -  u)dy  =   t 
dx     o  dx   o  o 


|_  (PU29)  +  P6.U  g  -  T,   .  a,  (f  )/n  PU2 


This  can  be  written  as 
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dx     dx  dx         Uo 
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ii  +  e  I  i-,  l(pui)  +Hidu 

dx     LpIT   dx        U  dx 


Again: 
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dx 

v  1  /n 
where    Q  =  «i  (jjq)' 

and      m  =  e 


pdx  = 


d(PU2)  +  H  dU 
-PU2      u 
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Hence,  it  is  seen  thats 

1  /   A  n+1 

(m6)i/n  ~   (m6)  =  m  n 
dx 


TTH+2 
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Integrating  yields s 
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Inspection  will  also  show  that  Equation  (7b)  could  have  been  obtained 
directly  from  Equation  (7),  using  the  relationship  between  the  two  inte- 
grating factors,  and  between  6  and  9  „  Note  that  C.  =  2(H  +  2)„ 
Substituting  for  m  in  Equation  (7b)  yields? 


_  _  n+1       .     .  .   x   TIH+2  n+1  „,_   n+1 

(pUH+2e)—  =  n_+_l  avl/n  j      (£U    —  ^  +  (    H+2Q  — 

Jx       TTl/n  o  o    o 
o  U 
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This  can  be  further  simplified  to  either  of  two  forms; 

n+1  „   _     n+1      r         .        .  ,      ,-x     n+1 

Q-r—        ,    TTH+2 ,,--—-        n+1         1/n    r 
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where 
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e(_!)i/*  =  p-n-    u~n 
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n  =  (H  +  2) 
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..        x     n+1 

S-i-ia/     p~Undx  +  C 
■     n  ■'x  o-1 
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1/2 


(E-12) 


For  the   case  of  negligible  density  gradient  Equation   (E12)    reduces 

(■op) 

to  Equation  (2208)  in  Schlichtingv   '. 
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for  which 


V  =  820  ft/sec 


p  *  .25  lbm/ft 

R  =  06T  ft 
Hence,  the  estimated  pressure  gradient  over  the  top  of  the  boundary  layer 
at  the  throat  is 


*E  =   vi 


T500^§  =  U.33?^ 
dn      R  ftJ        in 

The  side  walls  are  approximately  l"  apart „   At  a  point  midway  between 

them,  the  pressure  gradient  is  zero.   Assuming  the  pressure  gradient  is 

linear, 

-+U.33 


|E.  =  _uo33  +  8.66z 
dz 


z  =  1.0" 


And 


p  =  p  -  U.33z  +  U.33z2 


w 


w 


.  -   r      **.33z2    *+.33z3  J 


=  °w  -  -Tl 


Hence;    d  -  p   =  .71  psi 
w   ^  * 

This  is  of  sufficient  magnitude  to  offer  a  possible  explanation  of  the 
less  than  1  psi  difference  between  average  mainstream  and  measured  wall 
pressure  at  the  throat. 

It  is  noted  that  greater  rates  of  curvature  of  the  boundary  layer  occur 
upstream  of  the  throat  at  relatively  negligible  velocities.   Downstream, 
where  velocities  continue  to  increase  the  boundary  layer  has  very  little 
curvature. 


Gl. 

APPENDIX  G 

DISCUSSION  OF  INTERFEROMETRY  APPLICATION  TO  TWO-DIMENSIONAL 
NOZZLE  FLOW 

Extensive  discussions  of  the  theory  behind  interferometry  techniques 
are  numerous.   References  (lU)  and  (36)  are  included  for  the  reader  who  may 
desire  additional  detail. 

Basically,  an  interferometer  consists  of  a  monochromatic  source  of 
light  which  is  focused  into  a  beam  of  nearly  parallel  rays .   The  beam  is 
then  split  into  two  beams  traveling  separate  paths  of  identical  distance 
after  which  the  beams  are  rejoined  and  focused  on  a  viewing,  or  photographic, 
plate.   An  interference  pattern  of  parallel  light  and  dark  bands,  as  seen  in 
Fig.  12  is  formed  when  the  two  beams  are  rejoined  due  to  the  fact  that  the 
various  rays  of  each  beam  travel  slightly  different  distances  depending  on 
the  geometry  of  the  mirror  arrangement. 

One  beam  is  passed  perpendicularly  through  the  test  section.   To  com- 
pensate for  the  change  in  optical  path  length  caused  by  the  glass  side  walls, 
"compensating"  glasses  of  exactly  the  same  thickness,  must  be  placed  in  the 
path  of  the  other  beam.   During  testing,  the  only  change  which  occurs  is  the 
density  of  the  vapor  in  the  test  section  which  influences  the  velocity  of 

light  passing  through  it  according  to  the  relation,  for  two-dimensional  flow: 

SX 
o 

-J-   =  nl  "  n2 
where  n  is  the  index  of  refraction  of  the  gas  in  the  test  section,  t    is  the 

width  of  test  section  through  which  light  passes,  X   is  the  wave  length  of 

light  being  used,  and  S  is  the  number  of  fringe  shifts  which  occur  at  the 

point  in  question  due  to  a  change  in  the  index  from  n  to  n  . 

The  index  of  refraction  is  related  to  the  density  by  the  relationship 

Gp  =  n  -  1 

where  G  is  the  Gladstone-Dale  constant,  fixed  for  any  gas. 


G2. 

The  technique  used  in  this  investigation  was  to  first  photograph  the 
interference  fringe  pattern  for  no  flow  and  then  take  another  photograph  under 
test  conditions.   The  fringes  of  each  were  then  numbered  from  a  common  point 
which  effectively  represented  stagnation  conditions .   For  each  distance  x 
it  was  then  a  simple,  though  time  consuming,  matter  to  determine  the  fringe 
shift  S  and  thence  to  calculate  the  density  ratio .   Photos  used  in  calcul- 
ating the  density  curves  in  Fig.  15(a),  (t>)  and  Fig.  16(a),  (b)  are  shown 
in  Figs .  11  and  12,  respectively,   A  graphical  illustration  of  the  determin- 
ation of  S  for  test  1  (Figo  15)  is  given  in  Fig.  13 »  Test  2  was  done  more 
precisely  by  measuring  the  distances  to  each  fringe  in  the  no-flow  picture 
as  measurements  in  the  flow  picture  proceeded,   since  it  was  detected  that 
slight  departure  from  perfectly  even  spacing  of  the  no-flow  fringes  occurs, 
probably  due  to  minute  differences  in  thickness  of  the  optical  glasses. 
(Note  differences  in  Tables  3  and  U„) 

Although  interferometric  methods  did  not  succeed  in  producing  measured 
effects  of  condensation,  as  mentioned  in  the  text,  the  results  obtained  for 
non-condensing  flow  does  offer  a  recommendation  for  possible  further  appl- 
cation  in  this  direction.   One  important  unfavorable  aspect  of  it,  however, 
is  the  extreme  amount  of  time  and  effort  required  to  reduce  the  photographs 
to  meaningful  data.   Each  of  the  two  tests  presented  here  require  alignment 
of  the  interferometer  with  the  nozzle  (a  process  which  may  take  as  long  as 
two  days,  after  which,  however,  it  is  permanent) ^  photography,  development 
of  negatives  and  prints  (another  day),  measurements  from  prints,  and  calcul- 
ations based  on  the  measurements.   If  pressure  measurements  are  taken,  coor- 
dinating the  location  of  pressure  and  density  measurements  also  requires 
painstaking  accuracy. 

Such  a  laborious  process  in  reducing  the  data  is  still  predicated  upon 
a  two-dimensional  nozzle  and  becomes  prohibitively  difficult  for  three- 
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dimensional  flow.   Reference  36  involves  the  analysis  of  one   three-dimen- 
sional nozzle  photograph. 


APPENDIX  H 


HI 


COMPUTATIONAL  PROGPAM 


Part  1:  Comparison  of  Equations  of  State 

Part  2:  Isentrope  and  Enthalpy  Calculations 

Part  3:  Normal  Shock  Calculations 

Part  h:  Nucleation  Contours  Calculations 

Part  5:  Condensation  Calculations 


PART   1  H2. 

*M38bl-3  3  74  »FMb  »DtdUb  »     i»l»lUti  i<,mKL     M    uur  F 

*  XEQ 

C       LOMPAKlbON  OF  tOUATlONo  oh  blATc 

KtAD  4,  IUiVU 
L       LALbULA I i ON  UolNb  DCJNLUiV-l   Wcdd  kudih  hwumTIuim  ur  sIhIl 

R=  .0820  / 

V=  Vo*44./16.08 

T=TO/] .8 

A0=2. 82392 

d0=0.0b20  74 

C0=0. 133670 

A=0. 133333 

6=0.00396667 

C=14814.8 

ALPHA=0.0U00890 

OAMM/  =0. UU3D7 

K8WK  =  k*T/v  +  (BU*R*r-Ao-CU/  (  T*T  )  )  /  (  V*V)  +  (  cs*k*  T-a  )  /  I  v**3  )  +  A*a 
1LPHA/ (  V**b  >     +     (C*(l.     +GAi'"»IA*(  V*V  )  )/  (T*T*V**3  )  )*tXPF  ( -bAili'iA/  (  V*V  )  ) 

PBWR  =  PdWR*14.696 

PRINT  It  PdWR 
C       CALLOLAlIuN  UMNG  HART  I  in  AND  nuj  LOuAIIuim  uF  S  I  a  1  c 

V  -     V0 

T  =  TO 

dV=.00749D 

ATwO=-8. 9273631 

dTwO= .OuD^b?t76 

CTWO  =  -1:>0. V73d / 

ATHRc=. 16907819 

£THRt=-.OJ007046l 7 

CTHRE=. 0831424 

AF0UR=-. 000211245^ 

dF  IVL  =  1.9b65593*10.**(-t>) 

TC=b47.b 

R=  .24381 

PMH  =  R*T  /  (  V-dV  )  +  I  ATrtO+bTwO*T+CTwU*tXPF  (  -:>«47b*T/  7  C  )  )  /  (  (  V-dV  )**2 
1)  +  (ATHKc.  +  dTHRt*T  +  C  I  riKt*tXPF  (  b  .47b*T  /  7C  )  )  /  (  (  V-dV  )  **3  )  +  AFuuk 
2/  (  ( V-dV)**4)  +  bF  I  Vt*T/ (  ( V-dV  )**b ) 

PRINT  2.  PMH 
C      CALCULATION  USING  FLANS,  mJuATIGN  OF  STATE 
C       FIND  LOCATION  OF  Ztku  Re  MM  I  NOcK 

R=  1.8644 

T=T0/1 .8 

J=l 

K=l 

PP=PBWR/14.696 

DELP=.  1*PP 
11  VP  =  R*T/PF+(.082b  +  .001265/*lJF)*]JOO./(  l.01*T)**3) 

VP=     VF/62.43 

IF     ( VP-V  )     b  ,1b  ,2b 
b  IF      (  1.0UOl*VP-V)     6,6,15 

6  IF      ( J-2  )      7  .8,6 

7  J  =  2 

8  IF      (K-2 )      10,9,9 

9  DELP  =     ,b*DELP 

10  pp=     f    j_[j£lp 


25 
26 
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28 
29 
30 

15 

1 
2 
3 
4 


GO  T 

IF  ( 

IF  ( 
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FORM 

FORM 

FORM 

FORM 

END 
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.9999*VP-V 
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15.26,26 


J-2)  30.29.29 

=.5*DELP 

PP+DELP 

0  11 

14.696*PP 

T  3.  PP 

AT  (lt>H  PkciiUKc  (  duR  )  =  » 

AT  (15H  PRtbiURHMH)  -, 

AT  (  15H  PRESSURE ( P  )   =  . 

AT  ( F5.  1  ,  F5.3) 


F  1 2  .  7  ) 


Fl. 


I   ) 


F  1 1 .  7  ) 


HU. 
PART  2 


MAIN  FOR  ISENTROPE  DECREASING  P 

DIMENSION  MARK(4)  »XOUT (200 )  ,  Y(l),  TOL(l),  YMIN(l) 
1   .  DY (  1  ) 

COMMON  PZERO  .HZERO   ,TZERO 

CALL  ENT  (P.T.H.4) 
100   READ4,PZER0»TZER0.DELTAP,MP 
4     FORMAT  (3F12.0.I4) 

X  =  0. 

CALL  ENT  ( PZERO»TZERO.Y( 1 ) ,3 ) 

HZERO=Y ( 1 ) 

TOL( 1 )=.0001 

YMINd  )  =10. 

H=DELTAP/3. 

XOUTd  )  =X 

DO  2  I=2«MP 
2     XOUT( I ) =XOUT ( 1-1 )+DELTAP 

MARK( 1 ) =1 

MARK(2 ) =MP 

MARK (4) =0 

CALL  RUNGE  (  1 , X , Y ♦ TOL , YM I N ,H , XOUT  ,MARK ) 

GO  TO  100 

END 


0022 
TOTAL     0022* 


*      LABEL 

C      DIFFEQ  FOR  ISENTROPE  CALCULATION 
SUBROUTINE  DIFFEQ  (N»X,Y,DY) 
DIMENSION  Y<  1  ) >DY(1 ) 
COMMON  PZERO 
P=PZERO-X 
IF(P) 1  ,2t2 

1  CALL  EXIT 

2  CALL  ENT  (P,T.Y(  1  )  ,2  ) 
CALL  FOR  RHO  (PtT.RHO) 
DY(1 )=-l./RHO 

RETURN 
END 

0013 

TOTAL     0013- 


H5. 

f      LABEL 

SUBROUTINE  PRINT  ( N  ♦  XOUT , Y , DY . J ) 
!      PRINT  FOR  ISENTROPE  CALCULATION 

DIMENSION  XOUT( 200) , Y( 1 ) »DY( 1 ) 

COMMON    P2ERO.HZERO    ,TZERO 

P=PZERO-XOUT( J) 

PRATIO=P/PZERO 

CALL  ENT  (P,T,Y( 1 ) .2 ) 

TRATIO=T/TZERO 

CALL  FORRHO  (P.T.RHO) 

CALL  FOR  RHO  ( PZ ERO , TZERO , RHOZER ) 

Z=P/(T*RH0*1.8644 ) 

RRATIO=RHO/RHOZER 

G=  (P-PSAVE)*(T+TSAVE) /(T-TSAVE)/(P  +  PSAVE) 

GAMMA=G/(G-1. ) 

GR=Z*GAMMA/(1.-(1.-Z )*GAMMA) 

VEL=0. 

W=2.*(HZERO-Y(l))*in33.*980.7 

IF(W)  6»6.7 
7     VEL=SORTF(W) 
6     H  =  Y(1  )*. 0242048+142. 60 

PSAVE=P 

TSAVE=T 

IF  (J-l)  2»2,3 

2  PRINT5»PRATIO»TRATIOfRRATIO.P  » T » GAMMA »GR »H » Z » VEL 
RETURN 

3  PRINTl  ,PRATIO.TRATIO.RRATIO»P.T»GAMMA.GRtH»Z.VEL 

5     FORMAT  (120H1        P/PO  T/TO     RHO/RHOO  P 

1  T        GAMMA  EFF        GAMMA  H  Z  VEL 

2  ♦/36X.84H  ATM        DEG  K 

3  CAL/GM  CM/SEC. / / 3F 1 2 . 4  ♦  F 1 2 . 3 » F 1 2 . 1 , 2F 1 2 . 4 ♦ F 1 2 . 2. 

4  .F12.4.F12.1 ) 

1     FORMAT  ( 3F12.4tF12.3.F12.1 t 2 F 1 2 . 4 , F 1 2 . 2 . F 1 2 . 4 , F 1 2 . 1 ) 
RETURN 
END 

0035 

TOTAL     0035* 


LABEL 

SUBROUTINE  FORRHO  (P»TtRHO) 

R  =  l  .8644 

V=R*T/P-(82.5+1.2657*P)/(.01*T)**3.333 

RHO=l./V 

RETURN 

FND 

0007 
TOTAL     0007* 


H6. 

C      P  IN  ATM  T  IN  DEG  K   H  IN  ATM-CC/GM 
SUBROUTINE  ENT  (P,T,H,D 
DIMENSION  HZEROI  130)  ,TH( 130) 

HOFPF(P.T) =4.333*(82.6*P+1.265*P*P)/(.01*T)**3.333 
GO  TO  ( 1»2»3»4) iL 

4  J  =  0 

7     J  =  J  +  1 
C       HZERO  IN  BTU/LB  MOLE   TH  IN  DEG  R 
READ5,HZERO( J) ,TH( J) 

5  FORMAT  (2F10.0) 
TH( J)=TH( J) /1.8 

HZERO ( J)=HZERO( J)*.52081 

IF(TH( J) )6»6*7 
C      NOW  HZERO  IS  IN  ATM-CC/GM  AND  T  IS  IN  DEG  K 
C      ZERO  H  AT  0  DEG    GASEOUS 

6  JEND=J-1 

1  RETURN 

3     DO  33  J=1.JEND 

IF(T-TH( J) ) 32.33.33 
33    CONTINUE 

GO  TO  100 
32    IR=J-1 

H  =  HZERO(  IB)  +  (HZERO(  IB+l)-HZERO(  IB)  ) /( TH(  IB+1 )-TH(  IB)  )*( T-TH(  IB)) 
1   -HOFPF(P.T) 

RETURN 

2  DO  23  J=1.JEND,5 
H1=HZER0( J)-HOFPF(P,TH( J) ) 
IFIH1-H)  23,25.24 
CONTINUE 
J  =  JEND 
JH=J-5 

IF(JH)  100,100,29 
DO  27  K=JH,J 

HA=HZERO(K )-HOFPF(P ,TH(K ) ) 
IF(HA-H)  27,28,28 
HB  =  HA 
GO  TO  100 
IB=K-1 

T=TH(IB)+(TH(IB+1)-TH(IB) )*(H-HB) / ( HA-HB ) 
RETURN 
T=TH( J) 
RETURN 
PRINT101 ,P,T,H 

FORMAT  (/42H  T  OR  H  IS  OUT  OF  THE  RANGE  OF  HZERO  INPUT 
1  ,/  12X,2H  P,12X,2H  T,12X,2H  H./3F14.2) 
CALL  EXIT 
END 


23 


24 


29 


27 


25 

100 
101 


0047 


TOTAL 


0047* 


HT. 


MAIN  TO  PUNCH  ENTHALPY  TABLES 

DIMENSION  HH(6) 

CALL  ENT(P,T,H.4) 

Tl=30. 

DO  2  L=l »4 

Tl=Tl+60. 

T2  =T1+10. 
PUNCH7.T2 
FORMAT  (F17.1) 
P  =  0. 

DO  2  J=l »25 
P=P+.5 
T  =  T1 

DO  6  K=l  »6 
T=T+10. 

CALL  ENT  (P,T,HH(K) ,3) 
HH(K)=HH(K)*. 0242048+142. 60 
H  TO  BE  PUNCHE  D  TO  AGREE  WITH  DIN  I  E 
PUNCH3,P .HH 
FORMAT  (F6.1  .6F11.2) 
CALL  EXIT 
END 


IC  CAL/GM    0  H  AT  0  SOLID 

0022 
TOTAL     0022* 


H8. 

PART  3 


LABEL 

SUBROUTINE  DIFFEO  (N.X»Y,DY) 

DIMENSION  Y( 2) »DY(2 ) 

COMMON  PZERO 

5=1. 

P  =  X 

GO  TO  (  1  »2  ) »N 

P=PZERO-X 

c  =  -l  . 

CALL  ENT  (P*T,Y( 1 ) »2 ) 

CALL  FOR  RHO  (P.T.RHO) 

DY(1 )=S/RHO 

DY(2)=1. 

RETURN 

FND 


0015 
TOTAL     0015* 


MAIN  FOR  SHOCK  CALCULATIONS 

CALL  ENT  (P,T.H,4) 

PRINTl 

FORMAT  (120H     TOLL  P01         T01  PI  fl 

1  P2         T2        P02        T02    P02/P01     P1/P02 

2  ,) 
CALL  START(l) 

END 

0008 

TOTAL     0008* 


LABEL 

SUBROUTINE  EQ1(V»REM) 

COMMON  PZERO,TZERO»HZERO,TOLL»P2»T2tPl.Tl»Hl 

CALL  FOR  RHO  (PI  .T1.R1  ) 

Ul=SORTF( 2.*(HZER0-H1 ) ) 

U2=U1*R1/V 

P2=P1+R1*U1*U1-V*U2*U2 

H2=HZERO-U2*U2*.5 

CALL  ENT(P2»T2.H2»2) 

CALL  FOR  RHO  (P2tT2»R2) 

REM=R2-V 

RETURN 

FNO 

0013 

TOTAL     0013* 
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*      LABEL 

SUBROUTINE  ENT  ( P  ,  T ♦ H  .  L ) 
C       P  IN  ATM  T  IN  DEC  K        H  IN  ATM-CC/GM 

DIMENSION  HZEROf 130 ) ,TH( 130 ) 

HOFPF(P.T)=4.333*(82.6*P+1.265*P*P)/(.0l*T)**3.333 

GO  TO  (  1  »2»3,4) ,L 

4  J  =  0 

7     J=J+1 
C       HZERO  IN  BTU/LB  MOLE   TH  IN  DEG  R 
READ5,HZERO( J) »TH( J) 

5  FORMAT  (2F1O.0) 
TH( J)=TH( J) /1.8 

HZERO ( J)=HZERO( J)*. 52081 

I F(TH( J) )6,6»7 
C       NOW  HZERO  IS  IN  ATM-CC/GM  AND  T  IS  IN  DEG  K 
C      ZERO  H  AT  0  DEG    GASEOUS 

6  JFND=J-1 

1  RETURN 

3     DO  33  J=1»JEND 

IF(T-TH( J) )  32.33.33 
33    CONTINUE 

GO  TO  100 
32    IB=J-1 

H  =  HZERO(  IB)  +  ( HZERO(  IB  +  1)-HZER0(IB))/(TH(IB+1)-TH(IB))*(T- 
1   -HOFPF(P.T) 

RETURN 

2  DO  23  JM.JEND.S 
Hl=HZERO( J)-HOFPF(P,TH( J ) ) 
IF(Hl-H)  23.25.24 

23  CONTINUE 
J=JFND 

24  JH=J-5 

IF(JH)  100,100»29 
29    DO  27  K=JH,J 

HA  =  HZERO(K  )  -HOFPF ( P , TH ( K )  ) 
IF(HA-H)  27,28.28 

27  HR=HA 

r,0  TO  100 

28  IR=K-1 
T=TH(IB)+(TH(IB+1)-TH(IB))*(H-HB)/ (HA-HB) 
RETURN 

25  T=TH(J) 
RETURN 

100  PRINT101 ,P ,T,H 

101  FORMAT  (/42H  T  OR  H  IS  OUT  OF  THE  RANGE  OF  HZERO  INPUT 
1  ./  12X.2H  P,12X,2H  T.12X.2H  H./3F14.2) 

CALL  EXIT 
END 


TH( IB) ) 


0048 


TOTAL 


0048* 


H10. 

SUBROUTINE  PRINT  ( N , XOUT , Y . DY ♦ J ) 
DIMENSION  Y( 2 ) ,DYI 2 ) ,XOUT( 1 )  ,R(3),REM(3) 
COMMON  PZERO,TZERO,HZERO,TOLL,P2,T2,Pl,Tl  .HI 
GO  TO  ( 3,4)  ,N 

3  Pl=PZFRO-XOUT( J) 

CALL  ENT  (PI  ,T1  ,Y(  1  )  ,2) 

H1=Y(1) 
C      FIND  FIRST  GUESS  AT  RH02 

CALL  FOR  RH0(P1,T1 »"(  1  )  ) 

CALL  FOR  RHO(PZERO,TZERO,RZERO) 

R ( 1  )=1.01*R(  1  ) 

R  (  2  )  =  R  (  1  ) 

DP=. l*(RZERO-R( 1 ) ) 

CALL  EOl  (R(  1  ) ,REM( 1 ) ) 
11     R (2 )  =  R( 2 )  +  DR 

CALL  EQ1  (R( 2) .REM( 2  )  ) 

IF(REM(  1  )*REM(2)  )  10,10,11 
10    R(3 )  =  (R( 1  )+R(2)  ) *.5 

CALL  EOl  (R( 3) ,REM(3  )        ) 

L  =  l 

IF  (REM( 3)*REM( 1  )  )  12,13,13 

12  L  =  2 

13  R(L )=R(3  ) 
IF((P(2)-R(1))*2./R(1 1-TOLL  )  14,14,10 

14  rONTINUF 
CALL    START(2) 

4  IF     (Y(l)-HZFRO)     5,6,6 

5  XB=Y(2) 
HB=Y( 1 ) 
RFTI1RN 

6  P02  =  Y( 2  )-( Y( 2 )-XB)*( Y( 1 1-HZERO)  / ( Y(  1  )-HB  ) 
CALL  ENT  (P02,T02,HZERO,2) 
POUT=P02/PZFRO 

POUTA=P1/P02 

PRINT7,TOLL,PZERO,TZERO,Pl,Tl , P2 , T2 »P02 , T02 »POUT » POUT A 

7  FORMAT  ( F10.6,4(F10.3,F10.1 ) ,2F10.5) 
CALL  START(l) 

END 

0038 

TOTAL     0038* 


H10. 

SUBROUTINE  PRINT  ( N , XOUT , Y , DY , J ) 
DIMENSION  Y( 2 ) » DY ( 2 ) ,XOUT( 1 )  ,R(3),REM(3) 
COMMON  PZERO,TZERO,HZERO,TOLL,P2,T2,Pl»Tl»Hl 
GO  TO  (3,4)  ,N 

3  Pl=PZFRO-XOUT( J) 

CALL  ENT  (PI  ,U  ,Y( 1  )  ,2) 

H1=Y(1 ) 
C      FIND  FIRST  GUESS  AT  RH02 

CALL  FOR  RH0(P1.T1 .R ( 1 ) ) 

CALL  FOR  RHO(PZERO,TZERO,RZERO) 

R  (  1  )=1.01*R ( 1 ) 

R  (  2  )     =  R  (  1  ) 

DR=. 1* (RZERO-R( 1 ) ) 

CALL  EOl  (R( 1 ) ,REM( 1  )  ) 
11     R (2 )  =  R( 2 )  +  DR 

CALL  EOl  (R( 2)  ,REM( 2  )  ) 

IF(REM( 1 )*REM(2) )  10,10.11 
10    R(3  )  =  (R(  1  )+R(2)  )*.5 

CALL  EOl  (R( 3)  ,REM(3  )        ) 

L  =  l 

IF  (REM( 3)*RFM( 1 ) )  12,13,13 

12  L  =  2 

13  R (L)=R( 3  ) 
IF((R(2)-R(1))*2./R(1 1-TOLL  )  14,14,10 

14  CONTINUE 
CALL  START(2) 

4  IF  (Y(l)-HZFRO)  5,6,6 

5  XB=Y(2 ) 
HB=Y( 1 ) 
RETURN 

6  P02=Y(2 )-( Y( 2)-XB)*( Y( 1 )-HZERO)/( Y( 1 )-HB  ) 
CALL  ENT  ( P02,T02,HZFRO,2) 
POUT=P02/PZFRO 

P0UTA=P1 /P02 

PRINT7,TOLL,PZERO,TZERO,Pl,Tl  ,P2  ,  T2 »P02 » T02 iPOUT  ,POUT A 

7  FORMAT  ( F10.6,4(F10,3,F10.1 ) ,2F10.5) 
CALL  START(l) 

FND 

0038 

TOTAL     0038* 


Hll. 


SUBROUTINE  START(N) 

DIMENSION  MARK(4).XOUT(200).Y(2)»TOL(2).YMIN(2)»DY(2) 

COMMON  PZERO.TZERO.HZERO.  TOLLtP2.T2 

H=.05 

DO  1  J=1.2 

MARK( J) =1 
1      YMIN(J)=20. 

GO  TO(  100*200)  ,N 
100   READ4,  PZERO.TZERO.P1 ,TOLL 
C      SET  CONSTANTS  FOR  RUNGE  INTEGRATION  DOWN  THE  ISENTROPE 

X  =  0. 

TOL  (  1  )=TOLl_ 

CALL  ENT  ( PZERO.TZERO.Yf 1  )  .?  ) 

HZFRO  =  Y(  1  ) 

MARK (4) =0 
4     FORMAT  (4F12.0) 

XOUT( 1 ) =PZFR0-P1 

CALL  RUNGE  (  1  ,  X  i  Y , TOL » YM I N ♦ H , XOUT . MARK ) 

CALL  EXIT 
200   X=P2 

CALL  ENT  (P2»T2  tY(  1  )  ,3) 

MARK (2 ) =(PZERO-P2)*2  0.+2  0. 

TOL(2 )=1. 

Y(2  )=X 

XOUT( 1 ) =P2 

MAR=MARK ( 2) 

DO  7  K=2»MAR 
7     XOUT (K ) =XOUT ( K-l )+.05 

CALL  RUNGE  < 2 . X , Y , TOL , YM I N , H , XOUT , MARK ) 

CALL  EXIT 

END 


0031 
TOTAL     0031* 


SUBROUTINE  FORRHO  (P,T,RHO) 

R=1.8644 

V=R*T/P-(8  2.5+1.26  57*P) /(.01*T)**3.333 

RHO=l./V 

RETURN 

FNH 

0006 

TOTAL     0006* 


H12. 

PART  k 

*  LABEL 

C      PROGRAM  FOR  RATE  CONTOURS 

C      READ  SECTION 

DIMENSION  R(5)»  PP(5)»  FP ( 3  )  ,  P(3) 
101   READ1,  RHOA,RHOB,SURFA,SURFB»WMOL 

1  FORMAT  (5F12.0) 

READ1,  TSTARTt TEND.DELTAT 

READ40  ,R 
40    FORMAT  (5E12.2) 

PRINT2  2  .  RHOA.RHOB»SURFA,SURFB»WMOL»R 
22  FORMAT  (14H1RATE  CONTOURS.// 

1  6H  RHOL=,F8.3.  2H  -,  F9.6.10H  (T  DEG  C)./ 

2  8H  SIGMA  =,F8.3»2H  -.F9.6.10H  (T  DEG  C),/ 

3  7H  WMOL  =»F6.1  .// 

4  24H  TABLE  CONTAINS  P  IN  ATM.// 

5  19X.6H  RATF=»E12.3,4E18.3,10H   P  INF        #/ 

6  9H  T  DEG  K  ./ 

7  5H     =  ,/  ) 
M=(TEND-TSTART)/DELTAT 
T=TSTART-DELTAT 

DO  20  1=1 ,M 

T=T+DFLTAT 
C      COMPUTATION 

RHOL=RHOA-RHOB*( T-2  73.16 ) 

SIGMA=  SURFA-SURFB#( T-273.16  ) 

CALL  PINFIN  ( T.PINF) 

IF  (PINF)  20,20,74 
74    Cl=1.01325E6*1.0E16/( 1.38*T)  *PINF 

C1=L0GF(C1 ) 

C1=L0GF(        SQRTF(WMOL*2.*SIGMA/ (3. 1416*6. 025E23) )/RHOL)+2.*Cl 

C2=2.*SIGMA*WMOL/<32.062*1.01325E6*RHOL*T) 

C2=C2*C2*4.188*SIGMA*1.0E16/(  1.3  8*T  ) 

DO  19  J=l ,5 

IF(R(J))  19,19.21 
21    A  =  LOGF(R( J  )  )-Cl 

B  =  -2. 

C  =  C2 
C      CHOP  FOR  P 

P(  3  )=0. 

FP(1)=C 
11    P(3)=P<  3 )  +  3. 

IF(P( 3)-55.  )  66,66,  19 
66    IF(FP(1)*(P( 3)*P(3)*( A+B*P(3 ))+C))2.10,ll 

2  P ( 1  )=P( 3)-3. 
P(  2  )=P( 3  ) 

6     P  (3  )  =  (P (  1  )+P( 2) ) *.5 

IF  (P( 2 )-P( 1 )-.005*P(2) )  10,10,9 

9  L  =  l 
FP(3)=P(3)*P(3)»(A+B*P(3) )+C 
IF  (FPU  )*FP(  3)  )  4,10,5 

4  L=2 

5  P(L  )=P( 3) 
GO  TO  6 

10  IF  (P(3)-70.)  18,18,19 
18    PP ( J)=EXPF (P( 3)  )*PINF 


H13. 


19  CONTINUE 

20  PRINT  33,  T,  PP  ,  PINF 

33    FORMAT  (  F 1 0 . 2 , 9X , 5E 1 8 . 3  ,  El  1  .  3 
GO  TO  101 

END 


TOTAL 


0060 
0060* 


PART  5 


Hlfc, 


MAIN  FOR  CONDENSATION 


100 


15 


24 


25 


D 

C 

1H 

2D 

3 

4 
I 

C 
R 
R 
R 
R 
P 
R 
R 
R 
R 
R 
R 
R 
R 
R 
R 
R 
R 
R 
F 
F 
P 

1  . 
F 

1 

2 

3F 

48 

5  = 
P 
T 
I 
X 
X 
X 

H 

M 

M 

M 

D 

T 
D 
T 
D 

Y 


IMENS 
OMMON 
FG,  S 
ADX, 
SSAT, 
,RTIM 
NTERP 
ALL  A 
EAD  1 
EAD  1 
EAD  1 
EAD  1 
EAD  1 
EAD  1 
EAD  1 
EAD  1 
EAD  1 
EAD  1 
EAD  1 
EAD  1 
FAD  1 
EAD  1 
FAD  1 
EAD  1 
EAD  3 
EAD  3 
ORMAT 
ORMAT 
RINT2 
SURFB 
ORMAT 
Fl 
GAMA 
12.2, 
H  ALP 
.F18. 
ZERO  = 
ZERO  = 
NTERP 
OUT(  1 
OUT  (2 
OUT(  3 
=  0.01 
ARM  1 
ARK  (2 
ARK(4 
0  24 
0L(  I  ) 
0  8  1 
0L(  I  ) 
0  25 
(  I  )=0 
YMINt 


ION  Y(10).TOL(10),YMIN(10)  ,  XOUTt  100)  »MARK( 5 ) 

PZERO.  TZEROtWMOL,  GAMMA,  RHOL, 
URFA,  SURFB,  XI,  ALPHA,  AMACH,  X,  AREA, 
P,  T,  RSTAR,  TD,  AMASS,  AVGR,  ASTAR, 

INTERP.Z  ,YDROP,YRATE  , GAMA , GAMB , Y  ,  NPAR T 
ES 
=  1 
REAC(X»AREA,DADX) 

PZERO 

TZERO 

WMOL 

Z 

GAMA 

GAMB 

RHOL 

HFG 

SURFA 

SURFB 

XI 

ALPHA 

RTIMES 

T0L(4) 

TOLL 

AMACH 

CNTAM 

RZERO 
(F  20.0) 
(El  1.3  ) 

PZERO,  TZERO,  WMOL,  Z , GAM A  ,  GAMB , RHOL , HFG » SURFA 

XI,  ALPHA,  RTIMES,  TOL ( 4  )  ,  TOLL,  AMACH,  CNTAM, 
(25H1  CONDENSATION  IN  NOZZLES,  ///9H  PZERO  =  , 
0.3,/9H  TZERO  =  ,  F10.2 ,/7HWMOL  =  ,  F12.2,/5H  Z  =  »F12.3,/7H 
=,F12.2,/9H  GAMB  =   ,F12.6,/8H  RHOL  =  ,F12.3,/7H  HFG  =  , 
/9H  SURFA  =  ,  F10.2./9H  SURFB  =  ,  F11.3./6H  XI 
HA  =,F14.5,/11H  RATE»10.**,F8.2,/9H  TOL ( 4 ) 
8,/9H  AMACHl  =,F10.2,/8H  CNTAM  =,E12 
2115. PZERO 
l.S*( TZERO+273. 16) 
=  2 

)  =  1. 

)=AMACH  #AMACH 
)  =2.»XOUT( 2 ) 

)  =  1 
1=30 

)=n 
1=5, Q 
=  TOLL 
1  =  1,3 
=  0.1 
I  =  1.9 


ATM 
DEG, 


GM/CC 
BTU/LBM 
DYNES/CM 
DYN/CM-C 


PTS/CM3 
CM 


,3,/8H 


=  ,F18 
RZERO 


RZERO 


,  F16.5,/ 
,8,/7H  TOLL 
=,E12.3,//) 


I  )=0. 00001 


H15. 


WMOL=WMOL/Z 

Y ( 1 )=25. 132  7*6  2.4*CNTAM*       1 544 .*TZERO/ ( WMOL*PZERO ) 

Y (2 ) =Y( 1 )*RZERO 

Y(3 )=Y(2) *RZER0/2. 

Y(5)=l. 

Y(6  )  =  1  . 

PZERO=PZERO/2115. 

TZER0=TZER0/1. 8-273.  16 

NPART=1 

X  =  0. 

N  =  Q 

CALL  RUNGE  (N.  X,  Y,  TOL .  YMIN,  H»  XOUT,  MARK) 

READ  77.  L 
77     FORMATH2) 

IF(L )  101. 101  ,100 
101    CALL  EXIT 

END 

0072 

TOTAL     0072* 

*      LABEL 

C      SUBROUTINE  TO  COMPUTE  NUCLEATION  RATE 

SUBROUTINE  R ATEC ( PR  ESS  ,  TEMP ♦  SIGMA.  SSAT,  RATE) 
DIMENSION  Y(10) 

COMMON  PZERO,  TZERO.WMOL,  GAMMA.  RHOL» 
1HFG.  SURFA,  SURFB.  XI.  ALPHA,  AMACH,  X.  AREA. 
2DADX.  P,  T,  RSTAR,  TD,  AMASS,  AVGR»  ASTAR,  SSAT 
3  , INTERP ,Z .YDROP .YRATE.GAMA.GAMB.Y.NPART.RTIMES 
A=RSTAR 

B=10,**16/ ( 1 ,380*TEMP) 
C=-4»188  7*SIGMA*A*A*B 

13  D=PRESS*1.01325E6*B 

14  E  =  WMOL/ (RHOL*6.025E23  ) 
EE=2.*SIGMA/( 3.1416*WM0L )*6.025E23 
AA=2.*LOGF (D)+  LOGF(E)  +  . 5*LOGF ( EE  )  +  C 
AA=AA+2. 303*16. 

IF(AA)  16,16,17 

16  RATE=0.0 
GO  TO  12 

17  CONTINUE 
IF(AA-85.0)  18,18, 1Q 

19     AA=85.0 

18  RATE=EXPF( AA ) 
12     RETURN 

END 

0025 

TOTAL     0025* 


Hl6. 


*      LARFL 

C      SUBROUTINE  TO  COMPUT 
SUBROUTINE  DIFFEO  (N 
DIMENSION  Y(  10)  ,DY(  1 
COMMON  PZERO.  TZERO, 
1HFG.  SURFA,  SURFB,  X 
2DADX,  P.  T,  RSTAR.  T 
3  SSAT,  INTERP.Z  .YDR 
55    TEMP=Y(6 )* ( TZERO+273 
R=.08478/WMCL 
GAMMA=GAMA-GAMB* ( TEM 
FAC=1./ ( l.+IGAMMA-l. 
GAMMA=GAMMA*FAC 
GO  TO  (  1 »2  )  ,NPART 

1  DO  3  1  =  1,9 
3     DY( I )=0. 

DY(5)=-Y( 5) /(2./GAMM 
DY(8)=R*TEMP*DY( 5 ) /Y 
DY(6 )=Y ( 6 ) *( GAMMA-1. 
RETURN 

2  CALL  PINFIN  (TEMP, PI 
PRESS=Y< 5 )*PZERO 
SSAT=PRESS/PINF 

IF  (SSAT-1.)  30 

30  RATE=0 
RSTAR=0 
DY(9)=0. 
DROPG=0 
GO  TO  32 

31  SIGMA=SURFA-SURFB*(T 
RSTAR=2.*SIGMA*WMOL 
IF  (Yd  1-25.1327  )  65 

65  AVGR=1 ,3*RSTAR 
GO  TO  67 

66  AVGR=SORTF (2.*Y( 3)/Y 

67  CALL  RATEC(PRESS,  T 
CALL  DROPGC(PRESS,  T 

6  DROPG=DROPG/(Y(7 )*SQ 
VEL=SQRTF  (-2.*Y(8)* 
DY(9)=1./VEL 

32  CALL  AREAC(X»  AREA, 
COMLAM= ( GAMMA-1 .0 )/G 
DY( 1 )=25.1327*  RATE* 
DY(2 )=Y ( 1 ) *DROPG+DY ( 
DY(3)=Y( 2)*DROPG+12. 
YDROP=RHOL/AMASS*Y( 3 
YRATE=RH0L/AMASS*4.1 
DY(4)=YDR0P+YRATE 

DO  444  1=1 »4 
444   DY( I ) =0. 

DY(5 ) =Y ( 5  )  *<  (  (CONLAM 
1  ( (GAMMA-1 . ) /GAMMA+1 

DY(6)=Y (6 )*( CONLAM*D 

DY(7) =Y(7 )*{-( 1  ,-Y(4 
1    .5/Y(6 )*DY(6)     ) 


E  DERIVATIVES 

,  X,  Y,  DY  ) 

0) 

WMOL»  GAMMA.  RHOL. 

I.  ALPHA,  AMACH,  X,  AREA, 

D,  AMASS,  AVGR,  ASTAR, 

OP.YRATE  ,GAMA,GAMB,Y,NPART 

.16) 

P-273.16) 
)*(  l.-Z  )  ) 


A+X* (GAMMA-1. ) /GAMMA) 

(5) 

)*DY( 5 ) / ( GAMMA*Y( 5) ) 


NF) 


30,  31 


EMP-273. 16) 

/(82.06  2*1.0132  5E6*RHOL*TEMP*LOGF(SSAT) ) 

,65.66 


(  1)  ) 

EMP,  SIGMA,  SSAT,  RATE) 
EMP,  SSAT,  AVGR,  DROPG) 
RTF(GAMMA*8.314E7/WM0L*TEMP) ) 
980. 7E6) 

DADX  ) 

AMMA*HFG*WMOL/ (15  44.*1.8*TEMP)*778, 

AREA 

1  )*RSTAR 

56637*  RATE*AREA*RSTAR*RSTAR 

)*DROPG 

888*RATE*AREA*RSTAR**3 


-l./(l.-Y(4)))*DY(4 )-DADX ) /( 1.  ( l.-YI 4) )* 
./(GAMMA*Y(7)*Y(7) )      ) ) ) 

Y(4)+(GAMMA-1.)/GAMMA*(1.-Y(4))*DY(5)/Y(5)) 
)  )*DY( 5 ) / (GAMMA*Y(7 )*Y(7)*Y(5))- 


■  IT-1 


DY(8 ) =R*TEMP*DY (5 )/Y( 5 ) 

RETURN 
END 


PINFIN  FOR  C02      0  To  240  DEG 

BELOW  93   PINF*C*tXPF  l-HFG.'R/ T  ) 

ABOVl  93  AN  ECJUATlUivl  (jF  Thl.  foR 

BETWttN  DATA  POINTS.      DATA  TAKtN  FROM 

CHEMISTRY.  I NTERNAT T UNAL  CRITICAL  TABLES.  AND  MATHtbON  GAi  OATA 


K 

«Ao  UotU  AND  PINF  AT  93  UitU  TO  G t T  C 
P=C*fcXPF  I  A/T  )  WAS  UiEU  TO  I  inT  tKPOLAT  E 
HAUDdOOI*.  OF  PHYSICa  AiMU 
RIT  ICAL 


11 
10 


su 

lc 

A  = 
B  = 
GC 
IF 
A  = 
B- 
GO 
IF 
A 
B 

GO 
1 
A  = 
B  = 
GO 
IF 
Ax 
B  = 
GO 
A" 
E  * 
H 
RE 
EN 


OUT 
T-9 
03  3 
263 
0  1 
T-l 
492 
430 
0  1 
T-l 
085 
373 
0  1 
(  T- 
835 
333 
0  1 
T-2 
028 
367 
0  1 
727 
70. 
'10 
RN 


INE  PINFIN  ( T.PINF) 
3.11,1,2 


0 


40.)  3»3,4 


60.)  5,5,6 
1 


195.  )  7,7 
5 


16.6  )  9,9,11 


0 


1 

.** ( A+B/T ) 


H18. 

*      LABEL 

C     SUBROUTINE  TO  PRINT 

SUBROUTINE  PRINTIN.  XOUT ,  YOUT.  DY,  J) 

DIMENSION  XOUT(l).  YOUTI7)   »Y(10)  »DY(1) 

COMMON  PZERO.  TZERO.WMOL.  GAMMA*  RHOL » 
1HFG.  SURFA,  SURFB»  XI.  ALPHA,  AMACH,  X,  AREA, 
2DADX,  P,  T,  RSTAR,  TD,  AMASS,  AVGR,  ASTAR, 
3  SSAT,  INTERP.Z  ,YDROP»YRATE  , GAMA  ,  GAMB , Y  .  NPAR T 

GAMMA=GAMA-GAMB*YOUT<6)*( TZER 0+273, 16) 

GAMMA  =  GAMMA/(  1,+  (GAMMA-1, )*(  1,-Z )  ) 

IF  (J-2)   21.22.32 

21  AMASS=ASTAR*PZERO*YOUT( 5 ) *2 1 1 5 . *SQRTF ( GAMMA* WMOL / ( 32.2*1544. 
1  *(TZERO  +  273.16)*YOUT(6)*1.8 )  ) 

AMASS= 32. 2*454. 6* AM ASS/ (30, 48*30, 48) 

PRINT10 

PSTAR=YOUT ( 5 ) 

TSTAR=YOUT (6) 

GSTAR=GAMMA 

RETURN 

22  NPART=2 

DO  26  1=1, 9 

26  Y( I )=YOUT( I ) 
DO  27  1=1,3 

27  Y(  I  )=Y(  I  )*AMASS 
Y(7)=AMACH 

A=PSTAR/(Y(5)*AMACH  )*SQRTF( Y( 6)*GSTAR/ ( TSTAR*GAMMA) ) 
INTERP=3 
ARFA=A*ASTAR 

CALL  AREAC  (X.AREA.DADX  ) 
INTERP=2 
X0UT(3) =X 
DO  25  1=3.63 
25    XOUTt 1+1 )=XOUT( I )+.2 
32     CALL  AREACtXOUTt J) .  AREA,  DADX  ) 
AREA=AREA/ASTAR 
TEMP=YOUT( 6 ) * ( TZ ERO+273 . 1 6 ) 
CALL  PINFIN  (TEMP,  PINF) 
PRESS=YOUT ( 5 )*PZERO 
SSAT=PRESS/PINF 

SIGMA=SURFA-SURFB*(TEMP-273. 16 ) 
IF(SSAT-1.)  69,69,67 

67  CALL   RATECIPRESS,  TEMP,  SIGMA,  SSAT.  RATE) 
IFIYOUT (1 1-25.1327)  69,69,68 

68  AVGR=  SQRTFI 2.*Y0UT( 3) /YOUT(  1  )  ) 

CALL  DROPGCfPRESS.  TEMP.  SSAT,  AVGR,  DROPG) 
DTD=TD-YOUT( 6)*(TZERO+273.16) 
PPSF=2115.*PZERO 
TDEGR=1.8*( TZ ERO+273. 16) 

T0TALN=Y0UT(1 )*WMOL  *PPSF/(25.1327*62.4*AMASS*1544.*TDEGR) 
69    PRINT11,X0UT(J) ,AREA,Y0UT(5 )  .YOUT (6)»SSAT.Y0UT(7)  ,Y0UT(4)  • 

1  YDROP.YRATE.DTD.AVGR.RSTAR.RATE.TOTALN 
RETURN 

10    FORMAT  (120H       X     A/A*    P/PO    T/TO    S  SAT    MACH    MOIST  BY  GROWTH 
1Y  GROWTH,   NEW  PART    TD-T        AVGR       RSTAR       NRATE      TOTALN 

2  .//) 


H19. 


11 


FORMAT  <F7.3,F8.3,2F7.4,F8.3,F7.3»F8.5,2E10.3,F7.2»4E10.3) 
END 


0057 
TOTAL     0057* 


LABEL 

SUBROUTINE  DROPGC ( PRESS , TEMP ,  SSAT,  AVGR,  DROPG) 

SUBROUTINE  TO  COMPUTE  DROP  GROWTH 

COMMON  PZERO,  TZERO  ,  WMOL V , GAMMAV » RHOL » 


ALPHA, 
AMASS. 


1HFG,  SURFA,  SURFB,  XI 
2DADX-  P.  T,  RSTAR,  TD- 
3SSAT,  INTERP 

J=l 

P=PRESS*1.01325E6 

T=TEMP 

R=AVGR 

L  =  0 

D(  1  )=0 

D{ 2  )=0 

TD=T-1. 

TD=TD+1 . 

CALL  PINFIN1TD,  PINF) 

PD=1.0132  5E6*PINF*EXPF(R 

GASCON=8.314E7/WMOL 

CONLAM= ( GAMMA-1 . ) /GAMMA*HFG*7 

BETA  =  P/SQRTF(6.2  8  32*GASCOi\|*T) 

BETAD  =  PD/SQRTF(6.2832*GASCON*TD 

DROPG  =  X I /RHOL* (BETA-BE  TAD) 


AMACH.  X,  AREA, 
AVCR.  ASTAR, 


TAR*T/(TD*R)*LOGF(SSAT)  ) 

78./<1544.45*i.ci*T  IMP)* WMOL 


AA   = 
AAAr  .5* 
A  = 
GO  TO  (8,9! 

8  J  =  2 
IF (A)  10,10,9 

10  DROPG=0. 
PRINT11,  PD,P 

11  FORMAT ( 8E12, 3 ) 
RETURN 

9  1  =  1 

IF (A)  2,2,3 
!      1=2 

3  D( I )=TD 

D3=  D( 1 ) *D( 2  ) 
I F ( D3  )  4,1,4 

4  TD=.5*(D( 1  )+D( 2  )  ) 
L=L  +  1 

IFIL-20)  5,5,60 
60    IF(DROPG)  10,6,6 
RETURN 
END 


1.  -  BETAD*TD/(BETA*T ) 

CONLAM-1 . )*GArViA/ (GAMMA-1. )*(1. -BETAD/ BETA) 
AA  +  AAA  -  (l./XI  -  1.)     *ALPHA*( TD/T-1 . ) 
,  J 


T,A,AA,AAA   ,X,Y(4) 


H20. 


210 

220 

230 

250 

10 

20 

30 


35 

40 

50 


60 
70 


80 
90 
C 

C 

c 

1000 


1030 


1040 


LABEL 


FIRST  ORDER  DIFF.  EU.  ROUTINE  —  ADJUST;.  bTEP  oIZE 


RUNGEIN,  X,  Y 
Y  (  1  )  »  Y  M  I  N  (  1  )  » 
DY I  50)  ,  YA( 50  ) 


SUBROUTINE 

DIMENSION 

DIMENSION 

KBTWN  =  1 

KBIG  =  1 

KLOW  =  1 

NCOUNT  =  15 

J  =  MARK  (  1  ) 

MAX  =  MARKI2) 

L  =  MARK (4) 

IF  ( L)   210, 

LTEST  =  1 

GO  TO  ;  30 

LTEST  =  2 

NUM  =  L 

DO  250   I  =  1 

SUB!  I  )  =  TOL( 

IF  (MAX  -  J) 

RETURN 

A  =  XOUT ( J )  -  X 

B  =  ABSFI 1 .E-7  *  X) 

IF  ( A  +  B )   40  .  35» 

IF  ( A  -  B)   50.  50. 

J  =  J  +  1 

GO  TO  10 

CALL  PRINTIN,  XOUT, 

J  =  J  +  1  I 

10 

-  1.5*H) 


TOL,  YMIN,  H,  XOUT,  MARK) 
TOL(l),  oUB(50),  XOUT(l),  MARK(l) 
FA(50),  FB150),  FC(50),  YKEEPI50) 


210,  220 


N 
1/32.0 
20,  30 


30 


35 
60 


Y,  DY,  J) 


GO  TO 
IF  (A 
H  =  A 
L  '=  L 
GO  TO 
IF  (A 
H  =  . 


+  1 
1000 
-  3.*H 
5*A 


70,  70,  tiu 


90,  1000,  1000 


DO  RUNGE-KUTTE-MERSON  INTEGRATION 


XA  =  X  +  h/3. 

XB  =  X  +  .5*H 

CALL  DIFFEQ(N, 

X  =  X  +  H 

DO  1030 

YKEEPl I 

FA (  I  '  = 

YA(  I  :     = 

CALL  DIFFEQIN, 

DO  1040   I  =  1 


X,  Y,  DY) 

I  =  1  ,  N 
=  Y(  I  ) 
H*DY(  I  ) 
Y( I )  +  FA( I ) /3. 

XA,  YA,  DY  ) 
N 


YA( I )  =  Y( I )  +  FA( I ) /6. 
CALL  DIFFEQIN,  XA,  YA, 
DO  1C50   I  =  1,  N 
FBI  I  )  =  H*DY(  I  ) 


+  H*DY( I 
DY) 


/6, 


H21. 


1050 


1060 


1090 


1100 

1110 
1120 
1130 

1135 
1140 

1150 


1160 


1165 
1170 

1172 


1174 


1176 


1178 


1180 

1  190 
1200 

C 
C 

c 

100 


YA(I)  =  Y(I)  +  ,125*FA(I)  +  ,375*FB(I) 

CALL  DIFFEGHN.  XB.  YA,  DY  ) 

DO  1060   I  =  1 »  N 

FC  (  I  )  =  H*DY (  I  ) 

YA(I)  =  Y(I)  +  .5*FA(I)  -  1.5*Fd(I)  +  2.*FC(I) 

CALL  DIFFEQ(N.  X.  YA,  DY ) 

DO  1130   I  =  1 »  N 

Y(I)  =  Y(I)  +  FA(I)/6.  +  ,666666667*FC(  I  )  +  H*DY(I)/6. 

J  =  Y(  I  ) 

IF  (ABSF(U)  -  YMIN(I))   H30»  1090.  1090 

KLOW  =  2 

E  =  .2*ABSF(U  -  YA( I )  ) 

IF  (  E  -  ABSF( TOL ( I ) *U)  )   1110,  1100.  1100 

KBIG  =  2 

GO  TO  1130 

IF  (  E  -  ABSF( SUB( I )*U)  )   1130.  1120.  1120 

KBTWN  =  2 

CONTINUE 

GO  TO  ( 100,  1135)  ,  KLOW 

GO  TO  (  1180  ,  1140)  .  KBIG 

NCOUNT  =  NCOUNT  -  1 

IF  (NCOUNT)  1150,  1150.  1170 

PRINT  1160,   X  ,  H 

PRINT  1165.  (I,  Yd).  DY(I),   I  =  1,  N) 

RETURN 

FORMAT  (58H4STEP  SIZE  HALVED  15  TIMES  CONSECUTIVELY  SINCE  LAST  PR  I 
INT   /29H  PROGRAM  TERMINATED  AT    X  =   ,  E16.8,  8H,    H  =  ,  £16.8. 
2//3H   I  ,  13X,  4HYU),  16X,  5HDY1I),//) 

FORMAT  (13.  7X,  2(E16.8,  4X ) ) 

KBIG  =  1 

IF  (H  -  B)   1176.  1  172.  1172 

X  =  X  -  H 

H  =  ,5#H 

DO  1 174   I  =  1 t  N 

Y(  I  )  =  YKEEPt  I  ) 

KBTWN  =  1 

KLOW  =  1 

GO  TO  1000 

M  =  15  -  NCOUNT 

PRINT  1178.  M.  X,  H 

PRINT  1165.  (I.  Y(I),  DYU).  I  =  1.  N) 

RETURN 

FORMAT  (41H4STEP  SIZE  BECAME  TOO  SMALL  FuR  COMPUTER . /20H  IT  HAS  BE 
1EN  HALVED  ,  12,  21H  TIMEb  CONSECU T I VE L Y . / 29H  PROGRAM  TERMINATED  AT 
2    X  =   ,  E16.8,  8H,    H  -  ,  E16.6.//3H   I.  UX,  4HY(I).  16X. 
35HDY ( I ) .// ) 

NCOUNT  =  15 

GO  TO  (  1190  .  1200)  .  KBTWN 

H  =  2  .  *  H 

KBTWN  =  1 

KLOW  =  1 

CHECK  FOR  INTERMEDIATE  PRINT  OUT 
GO  TO  (10,  110 ) .  LTEST 


H21. 


1050 


1060 


1090 


1100 

1110 
1120 
1130 

1135 
114  0 

1150 


1160 


1165 
1170 

1172 


1174 


1176 


1178 


1180 

1190 
1200 

C 

c 
c 

100 


YA(I)  =  Y(I)  +  .125*FA(I)  +  .375*FB(I) 

CALL  DIFFEQ(N.  XB.  YA,  DY) 

DO  1060   I  =  1  .  N 

FC (  I  )  =  H*DY(  I  ) 

YA(I)  =  Y(I)  +  .5*FA(I)  -  1.5*Fd(I)  +  2.*FC(I) 

CALL  DIFFEQIN,  X,  YA,  DY  ) 

DO  1130   I  =  1  »  N 

Y(I)  =  Y(I)  +  FA(I)/6.  +  .666666667*FC( I )  +  H*DY(I)/6. 

J  =  Y(  I  ) 

IF  (ABSF(U)  -  YMIN(I))   1130.  1090.  1090 

KLOW  =  2 

E  =  .2*ABSF(U  -  YA( I )  ) 

IF  (  E  -  ABSFI TOL ( I ) *U)  )   1110.  1100.  1100 

KBIG  =  2 

GO  TO  1130 

IF  (  E  -  ABSFI SUB( I )*U)  )   1130.  1120.  1120 

KBTWN  =  2 

CONTINUE 

GO  TO  ( 100,  1  135 )  ,  KLOW 

GO  TO  ( 1180 .  1140) .  KBIG 

NCOUNT  =  NCOUNT  -  1 

IF  (NCOUNT)  1150.  1150.  1170 

PRINT  1160,   X,  H 

PRINT  1165,  (I,  Y(I),  DY(I),   I  =  1,  N) 

RETURN 

FORMAT  (58H4STEP  SIZE  HALVED  15  TIMES  CONSECUTIVELY  SINCE  LAST  PR  I 
INT   /29H  PROGRAM  TERMINATED  AT    X  =   ,  E16.8,  8H,    H  =  ,  E16.ti, 
2//3H   I  ,  13X,  4HYII),  16X,  5HDY1I),//) 

FORMAT  !IJ,  7X,  21E16.8,  4X ) ) 

KBIG  =  1 

IF  ( H  -  B) 

X  -     X  -  H 

H  =  ,5*H 

DO  1174   I 

Y  (  I  )  =  YKEEP (  I  ) 

KBTWN  =  1 

KLOW  =  1 

GO  TO  1000 

M  =  15  -  NCOUNT 

PRINT  1178,  M.  X,  H 

PRINT  1165,  (I,  Y(I),  DYU),  I  =  1.  N) 

RETURN 

FORMAT  (41H4STEP  SIZE  BECAME  TOO  SMALL  FuR  COMPUTER . /20H  IT  HAS  BE 
1EN  HALVED  ,  12,  21H  TIMEo  CONSECUT I VELY./29H  PROGRAM  TERMINATED  AT 
2    X  =   ,  E16.6,  8H,    H  ^  ,  E16.6,//3H   I,  l^X,  4HY(I),  16X, 
35HDY ( I ) ,// ) 

NCOUNT  =  15 

GO  TO  (  1190  ,  1200  )  ,  KBTWN 

H  =  2  .  *  H 

KBTWN  =  1 

KLOW  =  1 

CHECK  FOR  INTERMEDIATE  PRINT  OUT 


1176,  1172,  1172 


=  1 


N 


GO  TO  (10,  110 ) ,  LTEST 


H22. 


110 
120 


130 

140 
150 


IF  (L)   120.  120.  150 
PRINT  J  JO,   X,  H 

PRINT  140,  (I,  Y(I).  DY(I).   I  =  1.  N) 
L  =  NUM 
GO  TO  10 

FORMAT  (5H  X  =  ,  E16.8.  4X ,  4HH  =  ,  E16.8,  11X,  lHI,  13X,  4HY (  I  )  , 
1 16X,  5HCYI  I  )  /) 
FORMAT  (55X,  13.  7X,  2IE16.8,  4X  )  ) 
L  =  L  -  1 
GO  TO  10 
END 


1 
11 

10 
12 


13 
14 

lr- 

16 
17 


2 

20 

23 


LABEL 

INTERPOLATION     ROUTINE  FOR  COMPUTING  AREA 

DADX  IS  CONTINUOUS 

SUBROUTINE  ARE AC ( XX , AREA tDADX ) 

COMMON  PZERO,  TZERO,  WMOL,  GAMMA.  RHOL, 
1HFG,  SURFA,  SURFB.  XI,  ALPHA,  AMACH.XX,  AREA 
2DADX»P,T,RSTAR,TD,aMASo,hvGR»ASTAR, 
3SSAT,  I  JTLRP 

DIMENSION   DD(50),X(5  0),A(50),D(50).DA(50) 


AREAC   PART  1 


INTERP=1,  READ   AREA  X  AND  ASTAR 
SET  X,DADX  TO  CORRESPOND  TO  ASTAR 
ASTAR  MUST  BE  IN  AREA  DATA  LIST 
GO  TO  (  1,2.3. M  »  INTERP 
1=0 
1  =  1+1 

READ1G,  X ( I ) ,A ( I ) 
FORMAT ( 2F15.0) 
IF( A(  I  )  )  11,12,11 
READ10,  XX, ASTAR 
IEND=I-2 

IEND  =  LAST  VALUE  OF  I  USED  MINUS  1 
AREA=ASTAR 
DADX=0. 
FIND  IZERO 
1=0 
1  =  1+1 

IF( X (  I  )-XX)    13,14,14 
IZERO=I 
D(  IZERO)  =0. 
DO  15  I=IZERO,IEND 

D(I+1)  =  2.*(A(  1  +  1  )-A(  I  )  )/ ( X(  1+1 )-X<  I  )  )-D(  I 
DO   16  1=2, IZERO 
J=  I  ZERO- I +  1 

D(J)=2.*(A(J+1)-A(J))/(X(J+1)-X(J))-D(J+1) 
DO  1^     1=2, IEND 

DD(I)  =  (D(  I  )-D(  1-1  )  )  /  (X(  I  )-X (  1-1 )  ) 
PLUG   IA  INTO  DD 
RETURN 


H23. 


AREAC 


PART 


GIVEN  X    OUTPUT  IS  AREA  AND  DADX 
IFIXX-X (  I  )  )  20,23,22 
1  =  1-1 

IF( XX-X( I ) )   20,23,23 
I  B=  I 
IA=  1  +  1 
GO  TO  24 


H2k. 


22    1=1+1 

IF  (  XX-X (  I  )  )   25,23,22 
25    IA=I 

I  B  =  I  -  1 
i       FIND  AREA  AND  DADX 
ZL,  AREA  =  A(IB)  +  (  XX-X  (  lb)  )*<DD(  I  A) 

1   *(XX-X(  IB))*. 5  +  D(  IB)  ) 
302   DADX=(DD(IA)#(XX-X(IB))+U(Ia))/AREA 
RETURN 


AREAC   PART  3 


C 
C 
C 

C 

c 

C      GIVEN  AREA  FIND  X  AND  DADX 

C    IF    MACH  IS  LESS  THAN  1  LOOK  BEFORE  THROAT 

C       IF  MACH  IS  GREATER  THAN  1  LOOK  AFTER  THROAT 

3  I=IZERO 
IF(AMACH-1.)  30,31,31 

30  1=1-1 

IF  ( AREA-AI I ) )   32,32,30 
32    I  B=  I 

IA=I+1 
GO  TO  3^ 

31  1=1+1 

IF ( AREA-A ( I ) )   35,32,31 

35  I  A=  I 

I  B=  I  -1 
C    FIND  X  AND  DADX 
34     IF(  D(  IA)-D(  TB  )  )  36,37,3S 

36  SIGN=1. 

301   XX=  (-D(  IB)+SIGN*  SQRTF ( D (  I  B  )  *D <  I B )  -2 . *DD  (  I  A ) * ( A (  I B ) -AREA )  )  ) 
1  /DD( IA)  +  X< IB) 
IF(XX-X(  IB)  )   38,39,39 

38  SIGN=-1. 
GO  TO  301 

39  IF(XX-X(  IA)  )  302  ,302,38 

37  XX=( AREA-A( IB) )/D( IB)  +  X(IB) 
DADX  =  D(  IB  )  /AREA 

4  RETURN 
END 


H25. 


SUBROUTINE  AREAC  ( X . AREA  ,  DADX  ) 

COMMON  PZERO.TZERO.WMOL t GAMMA i RHOL . HFG » SURFA » SURFB t X  I ,AL»AM»XX,AREA 
1A,DADX»P,T,RSTAR,TD, AMASS , AVGR  ,  ASTAR » SSAT , INTERP 
GO  TO  (1  »2»3)  .INTERP 

1  CONTINUE 
ASTARali 
X=.55 
RETURN 

C      NOZZLE  TWO  B 

2  IF  (X-l  .45  )  10.11  ,11 

10  AREA=l.  +  .00889*(X-.55  ) 
DADX=.00889/AREA 
RETURN 

11  IF  (X-4. 85)  12,13,13 

12  AREA=1.008  +  .0218*(X-1.45  ) 
DADX=.0218/AREA 

RETURN 

13  AREA=1.082+.0133*(X-4.85) 
DADX=.0133/ARFA 

RETURN 

3  IF  (AREA-1.008)  30,30,31 

30  X= (AREA-1. )*112.5+.55 
DADX=.00889/ARFA 
RETURN 

31  X= (AREA-1.008) /. 0218  +  1. 45 
DADX=.0218/AREA 

RETURN 
END 

0028 

TOTAL     0028* 


II. 

APPENDIX  I 

ESTIMATION  OF  EXPERIMENTAL  ERROR 

TEMPERATURE  MEASUREMENTS 

Estimated  error  (difference  between  measured  and  actual  temperature) 
on  temperature  readings  except  where  otherwise  noted  is  +  0.5°,  -  1.5°F. 
This  was  determined  after  the  following  measurements  and  observations. 

1)  Temperature  of  various  portions  of  the  ice  bath  of  the  reference 
junction  were  found  to  vary  by  as  much  as  1°F  from  the  freezing  point, 
32.2°F,  despite  its  containment  in  an  insulated  thermos. 

2)  A  calibration  of  the  potentiometer  (Wheatstone  bridge  type)  used  for 
measurements,  compared  exactly  with  those  readings  of  the  standard  copper- 
constantan  conversion  tables  at  each  end  (32°F  and  212°F)  of  the  cali- 
bration and  showed  a  maximum  variation  of  +  1.0 °F  between. 

3)  The  time  constant  of  the  thermocouple  was  measured  with  the  aid  of  a 
Sanborn  Recorder  capable  of  recording  a  50  cps  signal  with  no  distortion. 
Subjecting  the  thermocouple  to  various  step  changes  in  temperature,  it 
was  found  that  for  still  water  t  =  .15  sec,  for  still  air,  t  =  l„i+  sec, 
and  for  response  to  the  variation  of  human  breath,  t  =  .6  sec.   The 
measured  time  constant  of  the  potentiometer  was  .9  sec  ±    .2  sec. 

Based  on  these  observations,  a  reference  temperature  of  32.2°F  was 
assumed   for  all  measurements  and  the  variation  in  ice  bath  temperature  al- 
lowed for  in  the  estimated  larger  possible  negative  error. 

Assuming  a  linear  temperature  transient  over  a  1  minute  period,  a  vari- 
ation from  1Q0°F  to  -  50°F  would  result  in  a  rate  change  of  -  U.0°F/sec. 
Hence,  for  the  non-condensing  curves,  this  effect  alone  could  contribute  an 
error  of  about  5°F,  in  the  thermocouple  (t  =  .6)  and  potentiometer  circuit 
(x  =  .9  sec).   Hence,  in  the  non-condensing  curves,  an  error  of  +  U.5°F, 
-  1.5°F  is  estimated.   In  measurements  involving  condensation,  care  was 
taken  to  let  temperature  readings  reach  a  relatively  constant  level  before 
data  (photos)  were  taken.   No  transient  effects  are  estimated  for  those 
cases . 


12. 
PRESSURE 

'KaaaBflBfi&cxgsgi  measurements 

Pressure  readings  taken  by  the  mercury  manometer  board  are  estimated 
accurate  within  t    .1"  Hg,  this  would  amount  to  t   1.25%  error  in  the  pressure 
ratios  near  station  18  on  the  figures  for  Nozzle  I. 

Minimum  pressure  readings  for  Nozzle  II  were  in  the  vicinity  of  28"  Hg, 
well  beyond  the  region  of  condensation,,   Hence,  pressure  ratios  in  this  vi- 
cinity are  estimated  accurate  to  within  „5%. 

In  the  region  of  condensation  (station  10)  pressure  levels  range  near 
80"  Hg.   Hence,  error  in  pressure  ratios  in  this  region  is  estimated  at 
±    .125%,  two  orders  of  magnitude  less  than  the  pressure  variation  caused  by 
condensation  at  this  point, 

Error  has  been  introduced  into  the  values  of  stagnation  pressure  wher- 
ever they  have  not  been  directly  measured  by  the  manometer  board  owing  to 
the  high  pressure  level, or  the  stagnation  tank  pressure  gauge,  as  occurred 
several  times  for  Nozzle  I.   General  observations  of  the  behavior  of  the 
pressure  ratio  at  the  throat  have  shown  that  although  its  value  remained 
remarkably  close  to  .535  for  all  pressure  levels,  it  had  a  slight  tendency 
to  decrease  with  increasing  pressure,   Based  on  these  observations  (and 
others  with  air,  where  average  p/p  measured  .523  and  the  same  tendency  was 
noted)  it  is  estimated  that  the  pressure  ratio  at  the  throat  could  have  been 
lowered  to  a  minimum  value  of  .525  at  the  highest  pressure  levels.,   Hence, 
calculated  stagnation  pressures  are  estimated  to  have  an  error  of  +  1% ,    -2% 

Measurements  in  Nozzle  II  used  the  stagnation  pressure  gauge  for  all 
determinations  of  stagnation  pressure  after  the  gauge  was  calibrated  and 
found  to  be  accurate  to  within  1/2  psi  at  all  pressures  within  the  range  of 
experimentation.   Consequently,  all  stagnation  pressure  readings  for  Nozzle 
II  are  estimated  accurate  to  within  1/2  psi  and  static  pressure  readings 
within  .1"  Hg, 


13. 

DENSITY 

Distance  between  optical  glass  side  walls  in  the  interferometric  meas- 
urements was  .956  inches.   At  the  wall  the  boundary  layer  was  about  .006" 
at  the  throat  and  somewhat  larger  downstream  (See  Figure  12).   Within  the 
boundary  layer  densities  are  lower  than  in  the  mainstream,  due  to  temperature 
recovery  effects .   This  will  cause  a  small  error  in  the  measured  interfero- 
metric fringe  shift  data  when  applied  to  determination  of  density  in  the  mid- 
stream.  An  upper  limit  on  error  in  density  measurements  was  estimated  by 
assuming  a  thickness  of  the  order  of  the  displacement  thickness  (also  given 
an  upper  limit  of  6*  =  1/3  6)  to  be  of  temperature  equal  to  the  stagnation 
temperature . 

Near  the  throat,  assuming  6*  =  .002" 

26*  =  .00V  I   =  ,956"  -   1.0" 

T  p 

stream    Q  6       Q 

— =  cO  »  =  .0 


Hence, 


stream 


p   -  p,  =  .2p  also,   p   -  p   -  „6p 

so      s  oss 


The  fringe  shift  due  to  difference  in  density  over  length  )i   is 


S=f  (p,-p2) 


Hence, 


S,   f    (ps  -  p{)    ,00,  C.2P,) 

S  "  I   (p  -  p  )  "  .96  (.6pJ  "  -1J/' 

OS  s 

Similar  estimates  further  downstream  indicated  that  maximum  error  in 
density  measurements  probably  do  not  exceed  .2%„ 
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